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Heat-transfer measurements were made for  filmwise 
condensation of steam on externally enhanced horizontal 
tubes under vacuum and at atmospheric pressure. Data were 
obtained for copper tubes with circular fins of rectangular, 
triangular, trapezoidal; and parabolic cross sections, for 
spiral fins of triangular cross section, for commercially 
available finned tubes and for wire-wrapped tubes. Four 
spirally finned tubes from each of Cu, Cu-Ni, Al, and stain- 
less steel and two tubes with fins of rectangular cross 
section from each of Cu-Ni and Al were manufactured and 
tested to investigate the effect -of thermal conductivity. 

Àmong spirally finned tubes, the optimum fin pitch was 
found to be 1.6 mm. The tubes with a parabolic fin shape 
showed the best performance with steam-side enhancements of 
4.1 and 6.2 under vacuum and at atmopsheric pressure, 
respectively. Enhancement ratios as high as 3.5 and 2.1l 
were obtained under vacuum and at atmospheric pressure, 
respectively, for the commercially available finned tubes. 
The heat-transfer performance decreased with decreasing tube 
metal thermal conductivity. 

For the wire-wrapped tubes, an optimum pitch to wire 
diameter ratio of about 5.1 was found, with steam-side 
enhancements of about 1.9 and 2.2 under vacuum and at atmos- 
pheric pressure, respectively. A recent theoretical anal- 
ysis of laminar film of  low-surface-tension fluids оп 
wire-wrapped tubes was modified to include the condensate 
retention of the tube due to the high surface tension of 
water. Agreement between this modified analysis and the 


experimental data was favorable. 
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I. INTRODUCTION 


À. BACKGROUND 

It is well known that the power required to operate a 
naval vessel at a given speed is proportional to its 
displacement. Therefore, a major effort is necessary to 
reduce the displacement in order to minimize the required 
power. One of the largest components of a naval vessel is 
the main condenser. In fact, present-day  condensers are 
equipped with smooth tubes, and therefore are large in size 
and weight. Increasing the performance or the effective- 
ness of the condenser can reduce the material and the 
construction cost and of course the weight. 

The effectiveness of the condenser is limited by the 
thermal resistances of the water side, the Steam side and 
through the tube wall. Generally, the thermal resistances 
of the water side and steam side are the most dominant. 
Reducing апу опе of these thermal resistances will 
contribute to an improved overall heat-transfer coefficient. 
Therefore, for a given heat duty, this corresponds to a 
smaller and lighter condenser. Improved  heat-transfer 
performance can be achieved by enhancement of the water side 
and/or the vapor side. Enhancement on the water side is 
possible with turbulence promoters, twisted-tape inserts, 
and deformation of the tube to produce a "roped" scheme, 
internal fins or ribs [1]. The main disadvantage of water- 
side enhancement is the requirement of increased power for 
pumping. Therefore,  vapor-side enhancement may promise 
better economic advantage, while the best advantage may be 
achieved by enhancing both sides based on a comprehensive 
analysis. The enhancement of the vapor-side can be achieved 
by using low-integral fins, roped tubes or fluted tubes or 


by applying coatings to promote dropwise condensation. 
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While externally finned tubes have been used since the 
1940s in order to enhance the vapor-side coefficient of 
tubes used in refrigeration systems, such tubes have not 
been used in steam condensers. The reason for this appears 
to be the common belief that externally finned tubes could 
not enhance steam condensation mainly owing to the large 
amount of condensate that floods between fins in the lower 
portion of the tube. Since the surface tension of water is 
four times greater than that of the refrigerants, a very 
significant proportion of the tube may trap water between 
fins, which could result in poor heat-transfer performance. 

The theoretical treatment of the steam condensation 
problem on horizontal finned tubes is very difficult due to 
the large number of controlling parameters, such as gravita- 
tional and surface tension forces, fin spacing, height, 
thickness and shape leading to three-dimensional flow of 
condensate. Due to the complexity of the problem, any theo- 
retical model requires simplifying assumptions which can 
lead to inadequate results. Therefore, a large pool of 
reliable data, systematically covering all of the relevant 
variables, is essential in order to test simplified theoret- 
ical models and/or to arrive at a satisfactory correlation. 

This thesis effort is a continuation of research being 
conducted at the Naval Postgraduate School (NPS) under a 
grant from the National Science Foundation. The basic test 
apparatus has been constructed by Krohn [2].  Graber [3] 
provided the instrumentation, and took preliminary data as 
the system experienced problems with non-condensing gases 
and partial dropwise condensation on copper tubes. Poole 
[4] made further improvements on the apparatus especially 
for leak tightness. He operated the apparatus both under 
vacuum and at atmospheric pressure, and tested a total of 
six finned tubes, with different fin spacing, as well as a 


smooth tube. Unfortunately, Poole had problems owing to the 
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occurrence of partial dropwise condensation. Using this 
system, Georgiadis [5] was able to obtain complete filmwise 
condensation on 26 copper tubes. The repeatability of data 
obtained by Georgiadis proved the accuracy of the test 
apparatus and associated instrumentation which was basically 
the same as that used by Poole [4] with some minor modifi- 
cations.  Georgiadis tested a total of 23 finned tubes with 
rectangular-section fins and three smooth tubes. He system- 
atically varied the fin spacing (0.5, 1.0, 1.5, 2.0, 4.0 and 
9.0 mm), fin thickness (0.5, 0.75, 1.0 and 1.5 mm) and fin 
height (1.0 and 2.0 mm). Table I shows the combinations of 
fin dimensions used for these tubes. Based on both vacuum 
and atmospheric runs, Georgiadis reported an optimum fin 
spacing of 1.5 mm and an optimum fin thickness of 0.75 to 
1.0 mm. Among the finned tubes with a fin height of 1.0 mm, 
the tube with a fin spacing of 1.5 mm and fin thickness of 
1.0 mm provided the best heat-transfer performance. This 
tube resulted in a steam-side enhancement (i.e., the ratio 
of steam-side coefficient for the finned tube to the value 
for the smooth tube at the same heat flux) of about 4 апа 
5.7 for vacuum and atmospheric pressure, respectively. He 
found that the heat-transfer performance was most sensitive 


to the fin spacing, while the effect of Fin thickness was 


relatively small. Further, he found that the performance 
increased with increasing fin height. However, he showed 
that the ratio Eo/Ar (i.e., the enhancement beyond the area 


enhancement) decreased with increasing fin height (for 
example, tube 6 with е = 1.0 mm gave Eo/Ar values of 2.13 
апа 3.01 for vacuum and atmospheric pressure, respectively, 
while tube 23 with e - 2.0 mm gave values of 1.69 and 2.25). 
It appears that the surface-tension induced thinning of the 
condensate film diminishes with increasing fin height. 
Continuing with this investigation, Flook [6] tested 19 
additional tubes (see Table I for details). These tubes 
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included two sets of four tubes with fin heights of 0.5 and 
1.5 mm, respectively. In addition, he studied the effect of 
fin shape using machined fins of triangular, trapezoidal, 
and "parabolic" fin shapes, (these tubes had a fin height of 
1.0 mm, a fin base thickness of 0.5 mm and a fin spacing of 
1.5 mm at the fin root).  Flook showed that the tube with 
parabolic fins (tube 38) outperformed the corresponding tube 
with rectangular-section fins (tube 17) by 10 and 15 percent 
under vacuum and at atmospheric pressure, respectively. As 
also pointed out by Flook, this tube did not have truly 
parabolic fins. Like previous researchers [7,8], Flook also 
pointed out that a fin shape, such as parabolic, that has a 
continuously decreasing curvature from fin tip to fin root 
provides considerable thinning of the condensate film, thus 
resulting in improved heat-transfer performance. 

Despite considerable achievements made by Georgiadis 
and Flook, the very complicated nature of the problem being 
Studied demands much more attention. This includes more 
testing to study the effect of fin shape, the effect of fin 
thermal conductivity, performance of commercially available 
tubes and the enhancement that can be achieved by wire- 


wrapping smooth tubes. 


В. OBJECTIVES 

The main objectives of this thesis are as follows: 

Take data on a number of tubes to check the repeat- 
ability with previous data [5,6], 

Take data on tubes with fins of different shapes (trian- 
gular, trapezoidal, parabolic, etc.), 

Take data on commercially available tubes, 

Take data on tubes with different thermal conductivity 
having rectangular, triangular, and spiral fin shapes, and 

Take data on wire-wrapped tubes with different spacing 


and wire diameter. 
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Develop a theory to predict the data for wire-wrapped 
tubes. 
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II. PREVIOUS INVESTIGATIONS OF FILM CONDENSATION 
ON EXTERNALLY ENHANCED HORIZONTAL TUBES — 


A. GENERAL OBSERVATIONS 

When vapor condenses on smooth horizontal tubes in a 
filmwise mode, the condensate flows down by gravity and a 
continuous film always exists around the tube. The latent 
heat released by the vapor will eventually be absorbed by 
the cooling liquid that flows through the tube. 

The condensate film resists this heat flow because of the 
low conductivity of the liquid. The resistance increases 
as the film thickness increases. At the top of the tube, 


the condensate film thickness is small and thereby the 


resistance is low and it increases with increasing 
distance around the perimeter of the tube. Since the 
thermal resistance of the condensate limits the heat- 


transfer performance of. the tube, to enhance heat transfer, 
it is necessary to reduce condensate film thickness. For 
horizontal tubes, thinning of the condensate may be achieved 
by uSing a finned, grooved or a fluted surface. 

In 1984, Yau et al. [9] measured the enhancement 
provided by copper finned tubes over smooth tubes for film- 
wise condensation of = steam. Similar experiments by 
Wanniarachchi et al. [10] also in 1984 confirmed that the 
observed enhancements were greater than could be explained 
by the increased surface area alone. This additional 
enhancement may be a result of the surface-tension forces 
which act to thin the condensate film. The effect of 
surface tension was first described by Gregorig [7] using a 
fluted surface. The surface tension induced a large pres- 
Sure gradient along the fin surface. This induced pressure 


gradient can be explained by using Figure 2.1. 
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Figure 2.1 Schematic of Condensate Profiie on 
Unflooded Fin. 


The pressure gradient due to the effect of the surface 
tension between a liquid and vapor is inversely proportional 
to the radius of curvature of the condensate surface. 

The pressure of the condensate at point A is higher than 
the vapor pressure because of the convex condensate surface 
at this point. The condensate surface at the valley is 
rather flat. This nearly infinite radius of curvature of 


the condensate surface results in no pressure difference 
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induced by the surface tension at this point. Therefore, 
the pressure at point B is almost the same as the vapor 


pressure. These pressures are given by: 


E om (2.1) 
РА ТР, | Y 
PBP EP (2:729 
B V 
where 

Py = vapor pressure, 

РА, Рв = liquid pressure at points A, B, and 

TA. FS ^ » radius of curvature of the condensate 


film at points A, and B. 


At point A, the radius of curvature is small, so the pres- 
sure at point A is higher than the pressure at point B (see 
equations (2.1) and (2.2)). Since, in reality, the radius 
of curvature changes along the condensate surface, between 
points A and B, the pressure within the condensate film 
varies along the height of the fin. The overall pressure 


difference between points A and B is given by equation 
ЖЕГІ” 
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where 


BEN pressure difference between points A and B. 


Since the radius of curvature of the condensate film at 
point A is very small, we can see from equation (2.3) that 
there is a large pressure difference between points A and B. 
This pressure difference causes the condensate to flow from 
point A to point B, thinning the condensate layer. On the 
other hand, the flow of condensate between the fins depends 
on the ratio of surface tension forces to gravity forces 
since the former acts to retain the condensate between the 
fins while the later acts to drain the condensate. As 
surface tension increases, the condensate tends to flood a 
lw eT теа Gf the eüUbe 1n which the condensate layer is 
thick and the thermal resistance increases, so a small heat 
transfer coefficient results. The flooded portion of the 
tube, as mentioned in section A, is defined by the retention 
angle, (y) (i.e., the angle from the bottom of the tube to 
the highest position of the tube where the interfin space is 
still full of condensate). The retention angle depends on 
the fin spacing, surface tension and gravity forces, and the 
fin shape. Therefore, on the one hand, using fins around a 
smooth tube increases the condensing area and thins the 
condensate film along the fin surface. However, these ђепе- 
ficial effects are offset by the flooding that occurs. 
Decreasing the retention angle increases the heat transfer 
performance. Therefore, any means reducing the retention 
angle is beneficial. As mentioned in section A, one way to 
decrease the retention angle is by attaching drainage strips 


at the lower part of the tube. 
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B. CONDENSATE RETENTION 

In 1946, the first measurements of condensate retention 
were made by Katz et al. [11]. These measurements were made 
under static conditions (i.e., no condensation taking place) 
uSing water, aniline, acetone, and carbon tetrachloride on 
ten different tubes with fin densities from 276 to 
984 fins/m, and fin heights from 1.2 to 5.7 mm. They meas- 
ured the retention angle by visual observation and by 
weighing the amount of retained liquid. Theoretical treat- 
ment of the problem using the measurement of surface tension 
by a capillary tube and by the pendant drop method was made 
to develop a formula to predict condensate retention as a 
function of condensate properties and the dimensions of the 
tube. Their result for the condensate retention angle, y is 


given by equation (2.4): 





v с r zl (2.4) 
siny eg | * (Dg7 - Do s 980 
4 
where 
gG = surface tension, 
pp = density of condensate, 
Б = acceleration of gravity, 
De = fin Diameter, 
5 = outside diameter of tube, and 
s = fin spacing. 


It was shown that condensate retention depends mainly on the 
ratio of surface tension to liquid density and on the fin 


spacing. 
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In 1981, Rudy and Webb [12] measured condensate reten- 
tion angles on three integral-fin tubes with three fin 
densities (748, 1024, 1378 fins/m). They used three 
different fluids (water, R-11l, and n-pentane) under both 
static and dynamic conditions. Their results showed that 
the retention angle increases with increasing surface 
tension to density ratio of the fluid. They also showed 
that the difference between static and dynamic retention 
angles was very small. For water, they reported that a 
Significant portion of the tube surface was flooded. 

In 1982, sRifert [13] reported equation (2.5) for the 
retention angle using a model of the capillary rise height 


of the fluid along a vertical plate. 


| 2 Бак (25) 
= os 1 - -e = 
Ņ = cos Pf бо А, | 
where 
P = wetted perimeter, 
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Later, in 1983, Rudy and Webb [14] developed an analyt- 


ical model to predict  condensate retention. They used two 


fin pitch, and 


H 


Profile area of the fin. 


finned sections, one in tubular form and the other by split- 
ting the tubular section and unrolling it into a vertical 
plate. They found that the vertical rise height of the 
condensate was the same for these two cases. Based on this 
observation, they modelled condensate retention on a flat 
plate to express the same on the finned tube. They made a 


simple force balance on the free body of condensate and 
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Schematic of Condensate Retention 


Figure 2.2 


on Finned Tubes. 


= 


developed an expression for the condensate retention angle 


as given by equation (2.6): 


у = cos"! | = ام‎ (2.6) 


Әсте 5 Ds 


where 
e = fin height. 


Both their analytical and test results showed that the 
retention angle increases with increasing surface tension- 
to-density ratio. Experimental results involving the use of 
water, e R-12, ammonia, and n-pentane were predicted to 
within 10 percent. 

Owen et al. [15] also recognized the necessity of 
including the effects of condensate retention in the heat- 
transfer models. The main simplifying assumption for their 
model was that the condensate retention angle was indepen- 
dent of condensation rate, so there is no difference between 
a static test and dynamic condensation. Therefore, they 
considered only a static analysis. A simple force balance 
between surface tension and gravitational forces resulted in 
an equation for the condensate retention angle as shown 


below: 


4 с (2.7) 


= сов“ ! | = 
у pr g s ОЕ s 


This equation is the same as equation (2.6), except that 
equation (2.7) is independent of fin thickness (t). A good 
agreement between this equation and the available data were 
reported by Rudy and Webb [12]. 
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In 1983, Honda et al. [16] performed experiments on 
finned tubes with and without porous drainage plates using 
R-113 and methanol as working fluids. They revealed from a 
photographic study that the static and dynamic profiles of 
the retained condensate were almost the same, and, by 
attaching a porous drainage plate, they demonstrated a 
Significant reduction in the retention angle. Considering a 
force balance between gravity and the surface tension force 
acting on the condensate, they made a theoretical analysis 


to predict condensate retention, leading to equation (2.8): 


4 0 сов0 | (2.8) 


= -1 - -------- 
ф = cos || pe B s Dg 


where 
0 - fin tip half-angle. 


They reported very good agreement between their theory, 
their own data and other available experimental data 
ТОТ 

Yau et al. [9] measured the condensate retention angle 
using water, ethylene glycol, and R-113 for finned tubes 
with and without drainage strips. They used an apparatus to 
simulate condensation on finned tubes. From their results, 
it appears that a drainage strip attached edgewise to the 
bottom of a tube has a significant effect on removing the 
condensate, so liquid retention is significantly reduced. 
They modified equation (2.8) in order to fit their experi- 
mental data, and developed the empirical relation given by 
equation (2.9): 


1.66 5 J 


ф = сов”! | - 
РЕ В $ ПЕ 
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§ = fin tip half angle. 


Continuing with their investigation on  condensate 
retention, Rudy and Webb [17], in 1985, modified their 
Previous model [14] for predicting condensate retention on 
horizontal tubes with fins of arbitrary shape. Experiments 
were made on four finned tubes with fin densities from 748 
to 1378 fins/m and one spine tube with a fin density of 1378 
fins/m. The fluids tested were R-1ll, n-pentane, and water. 
In addition, they tested a Thermoexcel-C tube with fin 
density of 1417 fins/m and R-1ll as the working fluid. As in 
the previous models, this model is based on the equal capil- 
lary rise height for a tubular section and another section 
that was made by splitting a tube section and unrolling it 
into a vertical flat surface. Equation (2.10) is recom- 


mended to predict the condensate retention angle 


2 И th) 


E | (92.510) 
ф = cos 1 - Dp ope [th +s) 6 А], 


where 
Ру = wetted perimeter of fin cross section, 
pu = fim base thickness, and 
Ap = profile area of fin over fin cross section. 


From equation (2.10), the retention angle increases for 
an increase of surface tension to density ratio of the 


liquid, fin density or for an decrease of tube diameter. 
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For the case of a horizontal tube with rectangularly-shaped 
fins, equation (2.10) reduces to equation (2.7). The exper- 
imental deviation from the predictive value of equation 
(2.10) was + 10 percent. 


C. THEORETICAL MODELS 

In 1948, Beatty and Katz [18] performed experiments with 
propane, n-butane, n-pentane, sulfur-dioxide, methyl chlo- 
ride, and Freon-12 condensing on Single finned tubes with 
fin densities from 422 to 630 fins/m to obtain the vapor- 
side heat-transfer coefficient. They used the Nusselt equa- 
tions for condensation on a horizontal tube and on a 
vertical surface, and considered the finned tube to be a 
combination of two parts, a horizontal plain tube and 
vertical fins. Thus, they expressed the average heat- 
transfer coefficient by a Nusselt-type equation based on an 
equivalent diameter. They modified the customary leading 
constant (0.728) found in the Nusselt equation to fit their 


experimental data and their correlation is given below: 


: " foe (or во а мај" E 
saco oao [ercer oo eme aan 
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(2.14) 
Aeff = Ar + n Ag gin) 
a ке: (2.16) 
where 
Aeff = effective area of finned tube, 
Ае - total surface area of finned tube, 
An = surface area of smooth tube, 
D, - equivalent tube diameter, 
hpk = average vapor-side heat-transfer coefficient, 
he, ="specific enthalpy of vaporization, 
ke - thermal conductivity of condensate, 
n - fin efficiency, 
AT - vapor-side temperature drop, 
Це - viscosity of condensate, and 
Py = density of vapor. 


The empirically determined leading constant (0.689) in 
equation (2.11) is only 5 percent less than the theoreti- 
cally derived constant (0.728) using Nusselt theory. But 
the average heat-transfer coefficient is greater than that 
predicted by Nusselt theory for a smooth tube since the 


equivalent diameter is smaller than the outside diameter of 
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the tube. They claimed a maximum error of + 7.2 percent and 
- 10.8 percent for the fluids they tested. However, they 
did not take into account surface-tension effects for thin- 
ning the condensate along the fin height. They also 
neglected condensate retention and assumed gravity to be 
entirely responsible for the flow of condensate. Over the 
decades following their work, many researchers have found 
this model to be quite adequate for low-surface-tension 
fluids and for tubes with moderate fin densities (i.e., for 
condensing fluid-fin density combinations to yield low 
retention angles). However, as the fin density or the 
surface tension increases, the model tends to overpredict 
the heat-transfer coefficient [12]. 

Some years later in 1971, analytical and experimental 
studies of condensation on horizontal tubes with trapezoi- 
dally shaped fins were performed by  Karkhu and Borovkov 
[19] for condensation assuming surface tension to create the 
dominant force. The analytical solutions were based on the 
following assumptions: 1) the thin condensate film repre- 
sents a laminar boundary layer; 2) surface tension causes a 
pressure gradient along the fin side; 3) gravitational and 
inertial terms in the equation of motion of the film along 
the side of the fins were small compared to surface tension 
terms and were neglected; 4) the motion of condensate in the 
trough area is laminar; 5) condensate drains by gravity into 
the trough; 6) no condensation takes place on the flooded 
portion of the tube; and 7) the fin temperature is constant 
along the height of the fin. Using Nusselt's basic assump- 
tions and the differential equation of condensate motion 
(assuming radial flow of condensate feeding into the 
interfin space) with appropriate boundary conditions, they 
were able to obtain the thickness of the condensate film in 
the interfin spacing (equation (2.17)). In order to calcu- 


late the temperature distribution along the fin height, they 
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assumed one-dimensional heat conduction. Using numerical 
methods to solve the resulting differential equations, they 


found expressions for the heat-transfer coefficient: 





Zp = 1.6 H 9*2 (1 - 0.35 H 79-3 m) (2.17) 
- _ С hfg 
a T (2.18) 
where 
2 A r Do (2.19) 
Pa ü р b+ cos 8 2 
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(2.20) 
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where 
b = half of fin tip width, 
e = fin height, 
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F. = effective condensate surface, 


5 
G = condensate flow rate, 
h = heat-transfer coefficient, 
Z = dimensionless depth of condensate 


between fins, 


dimensionless depth of condensate at fin base, and 


0 = fin semivertex angle. 


Experiments were performed with four different finned tubes 
to condense both steam and R-113 with slowly moving vapor 
and when the Weber number (equation (2.23)) is greater than 
10. Using the measured temperature at the fin root, they 
found equation (2.21) for the dimensionless depth of the 
condensate at the fin base within + 2 percent of the experi- 
mental data. Also, they solved the heat conduction equation 
over the fin to find the temperature distribution over the 
fin height (equation (2.19)). They found the vapor-side 
coefficients to be 50 to 100 percent greater than that for a 
smooth tube. Further, they reported that their predictions 
agreed to within +5 percent with their experimental data. 

In 1973, Edwards et al. [20] reported an analytical 
model for condensation on circumferential grooves on hori- 
zontal tubes that included the surface tension effect, 
gravity, viscous, capillary pressure, and condensate accel- 
eration during the flow around the tube. This model is 
based on the following assumptions: 1) the condensate pres- 
sure is uniform over any cross section; 2) the radius of 
curvature of the meniscus in the flow region ас the trough 
is constant; 3) the heat transfer and vapor friction on the 
meniscus are negligible; 4) the draining condensate from the 


fin side has zero velocity; 5) the grooves have small 
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height compared to the tube radius; and 6) the film has a 
contact angle to the fin tip. 

Using conservation of mass, with the overall heat- 
transfer coefficient as a function of local pressure differ- 
ence and making a simple force-momentum balance over an 
element of condensate film, they found a relationship for 


the local heat-transfer coefficient as given by equation 
125246): 
1/2 


2 2 E Ка өр (2.24) 
Өс + до | 


where 


1 1 
бот > [z (ln = + 0.11593] 2 сос Өр (2.25) 


i E kf Po (2.26) 
=© с. 


Ке = liquid thermal conductivity, 


КИ = Line thermal conductivity, 


w = groove width, 
Ос comeactmangle, and 
0, = groove half angle. 


They assumed no heat transfer through the flooded portion of 
the tube. Further, they assumed only circumferential flow 
of condensate, thus neglecting any flow along the fin 
surface in the radial direction. While this assumption 
could result in poor predictions, they did not provide a 


comparison of their theory with any experimental data. 
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In 1977, Zozulya et al. [21] modified their previous 
model [19] to find expressions for the rate of heat 
transfer. Using the differential equation of condensate 
motion at the trough, and the average temperature difference 
determined by numerical methods, they arrived at equation 
(2.27) which gives the dimensionless height of the conden- 


Sate at the interfin spacing: 


z = j anp o. (227) 


2 61/2 (к, р. ДТ) 3/* соз 5/ + Do 


Е; = (2.28) 
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b = thickness at top of fin, 

ср = thickness at fin base, and 

z = dimensionless condensate film thickness 


in the interfin spacing. 


They compared the results of equation (2.27) and available 
data of refrigerants (R-ll, R-12, and R-21) for condensation 


on finned tubes manufactured of copper, brass, and steel 


with rectangular and trapezoidal fin-shapes. Also, they 
compared experimental data for condensation of R-113 and 
steam on different finned tubes. Discrepancies within + 15 


percent were reported. 

In 1979, Webb [22] reported a procedure for the design 
and optimization of a fin surface for heat-transfer perform- 
ance. Equations (2.31), (2.32) and (2.33) were recommended 
to calculate the optimum profile given by equation (2.29) in 
order to maximize the heat-transfer coefficient given by 
equation (2.30): 
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where 


ћи - heat-transfer coefficient, 


P = projected area of convex surface, 
S = value of s at Q- On? 
0 = angular coordinate measured, 


from the crest of convex surface, and 


bee - Pics. 


According to the author, this model underpredicts the heat- 
transfer coefficient. The calculated augmentation ratio 
based on the projected surface area h,/ Ayu ranged from 3.4 
to 3.8 for tubes with length from 4 ft to 40 ft while his 
experiment showed values in the range from 4 to 8. 

In 1980, Rifert [23] studied condensation of stationary 
vapor on horizontal finned tubes enhanced by the effect of 
surface-tension forces that tend to pull the condensate to 


the fin root. In his analysis, he divided the tube into two 
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zones: a) the unflooded zone where the condensate film is 
thin, and b) the flooded zone where the condensate film is 
thick. He solved the two-dimensional heat-conduction 
problem for the wall by numerical methods for each zone and 
then determined the mean heat flux. In cases where conden- 
sate is retained in more than half of the tube perimeter, 
Rifert pointed out that a three-dimensional form of the 
heat-conduction equation must be used. Solutions to these 
equations by numerical methods revealed that, in most cases, 
the fin temperature is very nonuniform and it depends on the 
properties of the wall and the vapor and the heat flux. He 
stated that for the highly non-isothermal fin surface, the 
use of and average temperature drop from vapor to the outer 
wall temperature (AT) yields computed heat flux values that 
are very sensitive to AT. Since this is unacceptable, he 
recommended the use of the average heat flux for the 
computation. 

In 1981, Adamek [24] presented a method for the design 
of an optimum surface for condensing heat-transfer perform- 
ance. Similar to other researchers [7,8,22], Adamek recog- 
nized the importance of surface tension on the heat-transfer 
performance of finned surfaces. Since the dominant force on 
the crest is the surface tension, he neglected gravitational 
forces in this region. He derived equations for the conden- 
sate film thickness (equation (2.34)) and the wall surface 
profile by defining the curvature as a function of the 
distance along the surface (equation (2.35)). Using equa- 
tion (2.36) for the the curvature of the profiles of the 
Wave crest, and the necessity that the pressure within the 
condensate must decrease from wave crest to the trough, he 
defined a family of condensate surface profiles, whose 
curvatures are given by equation (2.38). A number of € 
values and their corresponding condensate surface profiles 


are shown in Figure 2.3. He found equations (2.37), (2.39), 


36 


and (2.40) for the film thickness, the average heat-transfer 
coefficient and condensate flowrate, respectively. The 
реввтесег E in equation (2.37) characterizes the aspect 
ratio of the fin cross section (ratio of the height to the 
thickness). As the aspect ratio increases, the parameter Е 
decreases. As shown by Adamek, € = - 0.5 represents the 
optimum surface for maximum values of the condensate flow- 


rate and the average heat-transfer coefficient. 


6 
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мћеге 
W(s) = wall profile, 
f(s) = film profile, 
к(5) = local curvature of the condensate surface, 
5 = length of path in liquid film, 
SD = length of convex surface over which 
the condensate flow, 
Š = film thickness, 
0 = rotation angle of normal to fin surface, 
o = angle from origin to Sn 
е - ratio of slenderness, and 
I - radius of curvature. 


In 1981, Shklover et al. [25] treated film condensation 


for finned tubes to investigate the effect of metal thermal 


conductivity on the heat-transfer performance. Stationary 
Steam was used as the condensing fluid. They showed that 
as the thermal conductivity decreases, the temperature 


difference through the film decreases and the temperature 
difference through the wall increases. For this reason, the 
finned tubes made of stainless steel or german silver have 
the same heat-transfer performance as the smooth tube, while 
brass and copper finned tubes outperformed the smooth tubes. 

Rudy and Webb [12] proposed a possible improvement to 
the Beatty and Katz [18] model by taking into consideration 


condensate retention. They applied equation (2.11) only for 
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Figure 2.3 Adamek [24] Condensate Surface Profiles. 


the unflooded part of the tube and recommended equation 


(2.41) for the average heat-transfer coefficient. 


= пи 
ћ = hax | = (2.41) 
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where hpk is computed using equation (2.11). But, since 
this is a gravity-based model and it neglects any heat 
transfer through the flooded portion of the tube, it under- 
predicts the average heat-transfer coefficient of condensing 
R-11 by 10 to 30 percent. They recognized that the surface 
tension effect must be taken into account. 

In 1982, Webb et al. [26] developed a new model which 
included surface tension effects. They modified the original 
Nusselt equation for a vertical plate so that surface 
tension causes the condensate to drain from the fin tip to 
the base and gravity causes the condensate to flow in the 
channel between the fins. Assuming surface tension as the 
dominant force along the fin side, they proposed equation 
(2.42) for the fin side coefficient he; | and the average 
heat-transfer coefficient for the entire tube is given by 
equation (2.43): 


1/% 1/4 


5 
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where h, was computed using the Nusselt [27] equation, and 


h = — hy + ng — Беін (2.43) 
р Р 
Тһів поде1 predicted the heat-transfer coefficient within 
+ 10 percent for R-12. 

Using equation (2.7) for the condensate retention angle 
they developed earlier, Owen et al. [15] modified the 
Beatty and Katz model to include the retention angle. They 
divided the tube into two parts: the unflooded portion, and 
the flooded portion with the condensation occuring on both 
the retained condensate and the fin tips. The equations 


necessary for this model are listed below: 
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where h, was computed using equation (2.11) with a leading 
coefficient of 0.725 instead of the value of 0.689, and 


ће = 


1 = 


1 
— + س‎ ° 
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where h, was computed using Nusselt [27] equation, and 


heff = 


КЕБЕ 


where 


heff ^ 


КеЕЕ = 


Kein = 


keff (2.46) 


e 


(1 - sN) kfin + sNkf (2.47) 


heat-transfer coefficient for the 


unflooded portion of the tube, 


heat-transfer coefficient for a 


plain tube, 


heat-transfer coefficient for the 


flooded portion of the tube, 


heat-transfer coefficient of the 


combined fin and retained condensate, 
effective thermal conductivity, 

fin thermal conductivity, 

number of fins per meter, and 


retention angle. 
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5 - fin spacing 


Owen et al. showed their model to predict the data for R-11, 
R-22, methyl chloride, n-pentane, sulfur dioxide, propane, 
and n-butane to within + 30 percent. However, as shown by 
Honda and Nozu [16], this model overpredicted the steam data 
by up to a factor of 2. 

In 1983, Rudy and Webb [29] developed a model based on 
surface-tension-induced linear pressure gradient along the 
fin height, thus assuming radial flow of the condensate 
feeding into the interfin space. Further, they assumed 
gravity-drained flow of condensate in the space between 
fins. The Nusselt equation for horizontal tubes was used 
for the tube area between fins, while the fin surface was 
treated by replacing the body-force term (i.e., "pg") in the 
Nusselt ‘equation by an equivalent expression based оп 
surface-tension-induced pressure gradient as developed by 
Webb et al. [26] and Rudy [28] earlier. Once again, they 
assumed no: heat transfer through the flooded portions and 


the resulting expression is given by equation (2.48): 


h = / 0.725 
bt Dg ug AT 
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where 
L - length of tube, 
A surface area of tube between fins, and 
Af, 7 fin ЗЕНОН Са 
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This expression provided an accuracy of better than 107 for 
condensation of R-11 on short, finely-spaced fins, but the 
accuracy dropped sharply with fins of increasing height and 
for larger fin spacing. This was, according to the authors, 
due to the assumed linear pressure gradient on the fin 
surface as this model is not valid when gravity forces 
become dominant (i.e., as e increases). Therefore, equation 
(2.48) is valid for fin densities from 1200 to 1400 fins/m, 
and fin heights of less than 1 mm. 

Continuing their research on film condensation, Honda et 
al. [16] did experiments on horizontal finned tubes by 
attaching a vertical drainage strip at the bottom of the 
tube to reduce condensate retention. Using R-113 and meth- 
anol as condensing fluids, they found vapor-side enhancement 
ratios (compared to the case without drainage strips) as 
high as 1.36 for R-113 and 2.08 for methanol. 

In 1985, Honda and Nozu [30] developed an analytical 
model for film condensation on horizontal low integral-fin 
tubes. They divided the tube into flooded and unflooded 
regions. This model is based on the following assumptions: 
1) the wall temperature is uniform; 2) the flow is laminar; 
3) the condensate film thickness is small; 4) the dominant 
flow on the fin is in the radial direction. Based on these 
assumptions, expressions for Nusselt number representing the 
flooded and unflooded regions were found. The average 


Nusselt number is given by equation (2.49): 


[О = Mugu nul = т 00: + Моағ пе(1 – Тре] — ФЕ) (2.49) 
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where 


pu - Fin Efficiency, 
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пр = fin efficiency for the flooded region, 

Мид = average Nusselt number, 

Мачы = average Nusselt number for the unflooded region, 
Nug¢ = average Nusselt number for the flooded region, 

Фф = angular coordinate, 


Ти. = dimensionless temperature for the unflooded 


region, and 


Tyg = dimensionless temperature for the flooded region. 


Comparison of the results of this model with the available 
experimental data showed agreement to within + 20 percent 
for 11 fluids and 22 finned tubes. However, their model 
overpredicted steam data by up to 40 percent. 

In 1985, Rudy and Webb [31] modified their previous 
models taking into account surface-tension effects on film 
drainage and condensate retention. They treated the conden- 
sation problem considering two major regions: unflooded and 
flooded regions. They further divided the unflooded region 
into finned area and the interfin area. They computed the 
average heat-transfer coefficient for the entire tube as 


given by equation (2.29). 





A Ағ aW y 
hrw = h no F [hr gÉ * n hen x JF T | he = (2.50) 


In order to compute the heat-transfer coefficient for the 
finned area (h¢;,), they used an expression developed by 
Adamek (equation 2.39)). One of the profiles they used for 
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the trapezoidal-shaped fin is shown in Figure 2.4. They 
used equation (2.38) with 0, - 85 degrees and for each fin 
geometry an iterative procedure to establish the E value for 
each profile to correct the fin thickness at the fin base 
given by equation (2.51). 


t = tp- tt + 2 5(Sm) (225) 


Further, since they assumed that the length of the convex 
surface is from fin tip to fin base, they corrected for the 


film thickness equation (2.52) resulting from the additional 


condensation at the fin tip: 


_ Би = За + 92-6 (2554) 
where ) 
K 
ы = ТЕ (2.53) 
r 


where h, is the heat-transfer coefficient for the interfin 
area in the flooded region. In order to calculate the 
interfin area (= they used the Nusselt equation with an 
iteration to account for the additional condensate drainage 
from the fins. Finally, to compute the heat-transfer coef- 
ficient in the flooded region they used the following 


equations: 


8 с ар2 / dbi (2.54) 
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пе = Ж ш (2.56) 


where dy; is the heat flux if the fin thickness were zero. 
Using a numerical technique, they computed 42 based on 
two-dimensional conduction through the fins and the conden- 


sate film. 
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Figure 2.4 Fin Geometry for the Webb et al. Model [31]. 
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D. FILMWISE CONDENSATION ON WIRE-WRAPPED TUBES 

Similar to fin tubes, surface-tension effects can be 
beneficial for filmwise condensation on wire-wrapped tubes. 
However, the surface-tension effect on wire-wrapped tubes is 
different than that on finned tubes mainly because of the 
very little heat transfer through the wires compared to fins 
that would transfer the majority of the heat. 

As can be seen from Figure 2.5, the condensate surface 
in the space between wires on the wire-wrapped tube is 
rather flat and the pressure difference between the conden- 
sate and the surrounding vapor is zero. However, due to 
the existence of a concave condensate surface at the point 
of contact between the wire and the condensate surface, a 
reasonable pressure difference within the condensate will 
exist from the inter-wire space to the immediate vicinity of 
the wire. In fact, the pressure at point A is the smallest, 
thus resulting in a condensate flow toward the wire, at 
which point the condensate would rapidly flow around the 
tube and eventually leave the tube at the lowest point. 
This axial flow of condensate toward the wires generally 
results in a smaller film thickness than on a plain tube. 
On the other hand, the film is generally quite thick in the 
region between the wire and the tube, resulting in negli- 
gible heat transfer through this region. As can be seen from 
Figure 2.5, the radius of curvature of the condensate 
surface near the wire strictly depends on the wire diameter. 
Further, the extent of thinning of the condensate film 
depends on the wire spacing. Based on the above-mentioned 
observations, an optimum combination of wire diameter and 
pitch must exist to yield the best heat-transfer 
performance. 

In 1985, analytical and experimental studies of conden- 
sation on horizontal wire-wrapped tubes were performed by 


Fujii et al. [32]. Their model is based on the following 
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assumptions: 1) the liquid film is very thin compared to 
the diameter of the tube; 2) the condensate film flow is 
laminar; 3) the inertia forces in the momentum equation of 
the condensate are negligible; 4) the shear stress at the 
vapor-liquid interface is negligible; 5) the properties of 
the condensate are constant; and 6) the temperature of the 
cooling surface is uniform. They also assumed that heat 
transfer occurred only through the thin film between wires 
and they neglected any heat transfer through the wires. 
They developed semi-theoretical equations to predict a heat- 


transfer coefficient enhancement ratio as shown below: 


Nuw 5 F2(A) 
LEV. Wo p ove 2 57 
Nus (s + Dw) F2(0) | | 
where 
1 (7 
Fo(A) = = 0] Fi1C6A) de (2790188) 
апа 


(tan( 9/2)) ^/3 (віпф)1/3 


күлле а а (2.59) 
| је а *(sing) ж 
0 


where 


4 g Dg 


= 2.60 
А EE ( ) 


2 д |3/2 
= =2 
А Dy (2.61) 
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es noticed that for A > 15 


& (1 +A) yh 
Р1(Ф,А) = Р](0,А) = лыы. | (NOD 
Therefore, in this situation with A » 15, 
4 (1 +A) Û A 
F2(A) | w > (СЕН) 


1/% 


"n (2.64) 


and equation (2.57) simplifies tol 


ao = ($) 


| 4 + A) | 1/* | 
Nu g 5 | (1 4 Wi . (2.65) 
Nus 


(s + Dy) 3 


They selected a value of 0.03 for the coefficient C in equa- 
meron (2.61) in order to fit their experimentel data. 
Experiments were performed with three wire-wrapped tubes to 
condense ethanol and R-11 vapor. Enhancement ratios of 3.7 
and 3.3 for R-11 and ethanol, respectively, were reported 
(see Figure 2.6). The maximum enhancement ratio occurred at 
a p/Dw value of 2. 


nn the original  Fuj et al. paper 
(2.65) contained an error y meee sama the erm (4/3) ji? 
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Figure 2.5 Condensate Film Profile on Wire-Wrapped Tubes. 
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Есһапо1 





Figure 2.6 Comparison Between Experimental Data 
and Semi-Theoretical Model of Fujii et al. [32]. 


III. DESCRIPTION OF TEST APPARATUS 


А. TEST APPARATUS 

The test apparatus used for this investigation was 
essentialy the same as used by Georgiadis and Flook [5,6]. 
A schematic of this apparatus is shown in Figure 3.1. Steam 
was generated using distilled water in a 304.8 mm (12 in.) 
Pyrex glass section which was fitted with ten 4000-Watt, 
440-Volt Watlow immersion heaters. The steam from the 
boiler flowed upward and passed through a 304.8 mm (12 in.) 
to 152.4 mm (6 in.) reducing section to a 2.44 m (8 ft.) 
long section of Pyrex glass piping. The steam flowed 
through a 180-degree bend and entered a 1.52-m-long section 
before finally entering the stainless-steel test section, 
which is shown in Figure 3.2. The test tube was mounted 
horizontally in the test Section. A portion of the steam 
condensed on the test tube, while the excess steam travelled 
downward and condensed in the auxiliary condenser. The 
condensate drained back to the boiler by gravity, completing 
the closed-loop operation of the system. | 

The exit side of the test tube was provided with a 
mixing chamber for accurate measurement of the outlet 
temperature of the coolant. A view port was provided in 
the test section to allow visual observation of the conden- 
Sation mode to ensure complete filmwise condensation during 
data collection. The auxiliary condenser consisted of two 
9.5 mm (3/8 in.) diameter water-cooled copper tubes heli- 
cally coiled to a height of 457 mm (18 in.). The auxiliary 
condenser was cooled by a continuous supply of tap water 
through a flow meter. A throttle valve was provided to 
control the flow rate through the auxiliary condenser, thus 
keeping the system at the desired internal pressure. For 


example, when the flowrate through the test tube was 
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decreased, the flowrate through the auxiliary condenser had 
to be increased. Filtered tap water was collected in a 
large sump with a capacity of about 0.4 cubic meters (Figure 
3.3), and was used to cool the test tube. Two centrifugal 
pumps connected in series took the water from the sump and 
passed it through a flow meter into the test tube. A valve 
on the discharge side of the second pump, and before the 
flow meter, allowed the velocity of water flowing through 
the test tube to be varied from 0 to 4.4 m/s (14.4 ft/sec). 
À vacuum pump was operated continuously during the oper- 
ation of the apparatus to remove non-condensing gases from 
the test section. The system used to remove non-condensing 
gases is shown in Figure 3.3. It was unavoidable that the 
vacuum pump mainly drew steam with trace amounts of air 
(non-condensing gases). To minimize the contamination of 
the pump by the steam, another condenser was provided to 
condense as much steam as possible. This condenser was 
cooled with the filtered tap water before it entered the 
large sump. The condensate from this steam collected in a 


Plexiglas cylinder to be drained later. 


B. INSTRUMENTATION 

The electrical power input to the boiler immersion 
heaters was controlled by a panel-mounted potentiometer. In 
order to compute the input power to the boiler, a root-mean 
converter with an input voltage of 440 VAC generated a 
signal which was fed to the data acquisition system. А 
more-detailed description of the boiler power supply is 
provided by Poole [4]. The temperatures of the steam, 
condensate and the ambient surroundings were measured using 
calibrated copper-constantan thermocouples made of 
0.25-mm-diameter wires. Two of them were used for the steam 
temperature, one for the condensate return and one for the 
ambient temperature. These thermocouples had an accuracy 


within + 0.1 К. when compared against a platinum-resistance 
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thermometer. Since the temperature rise of the coolant 
through the test tube is the most critical measurement, 
considerable attention was paid to obtaining the highest 
possible accuracy. For this purpose two independent temper- 
ature measurement techniques were used: a Hewlett-Packard 
(HP) 2804A quartz thermometer with two probes having an 
accuracy of + 0.02 K, along with a 10-junction, series- 
connected copper-constantan thermopile with a resolution of 
0.03 K. 

For most of the data collected, the quartz thermometer 
and the thermopile agreed to within + 0.03 K and when the 
difference was more than + 0.05 K, the data set was disre- 
garded and a repeat set was made. The cooling water flow 
rate was measured using a calibrated rotameter and the value 
was fed manually to the computer. Another rotameter was 
Provided to measure the cooling water flow rate through the 
auxiliary condenser. 

À pressure tap located about 50 mm above the test tube 
was connected to a U-tube, mercury-in-glass manometer gradu- 
ated in millimeters to measure the absolute pressure of the 
system. At the beginning and at the end of each test run, 
an accurate pressure reading was made and entered into the 
computer. The measured system pressure and the saturation 
pressure corresponding to the measured steam temperature 


were used to compute the concentration of any air that might 


have been present. For this purpose, a Gibbs-Dalton-type 
relationship was used. The computed non-condensing gas 
concentration was found to be within - 1.5 to QO percent. 


Such a value revealed that major air leaks did not take 


place following the last vacuum test on the apparatus. 


C. VACUUM INTEGRITY 

Vacuum tightness for any condensation system, especially 
at low pressures similar to marine-vechicle condensers which 
operate at ап absolute pressure of 2 inHg, is very 


important. 
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The reason for this is because any small amount of air or 
other non-condensing gas present with the condensing vapor 
tends to accumulate at the liquid-vapor interface. When 
this phenomenon takes place, an added thermal resistance 
occurs at the interface, which will degrade the heat- 
transfer performance considerably. Therefore, in order to be 
able to collect consistent and reliable data, exteme care 
was taken to ensure a leak-tight apparatus. In fact, during 
the early stages of this investigation, a major leak was 
found through the screws at the test section. After this 
was fixed, a vacuum test was carried out routinely about 
once a month. A leak rate that corresponds to a pressure 
rise of about 2 mmHg in 24 hours nearly at the operating 
pressure was found. Also, since the vacuum pump was oper- 
ated continuously during the experiment, any accumulation of 
non-condensing gases within the apparatus was effectively 


eliminated. 


D. DATA ACQUISITION SYSTEM 

An HP-9826A computer was used to control an HP 3497A 
Data Acquisition System to monitor the system temperatures 
and boiler input power (using the converter signal). Raw 
data were processed immediately and stored on diskette for 
reprocessing at a later time. After all the runs were 
collected, the data were reprocessed using a new Sieder-Tate 
Coefficient found by the modified Wilson method. 


ECOWNEUBES TESTED 

For this thesis effort, a total of twenty six tubes were 
manufactured. Table I lists all the finned tubes tested and 
their dimensions, including two tubes tested also by 
Georgiadis [5] and Flook [6] and four tubes tested by Flook. 
Twenty one of them were made from copper, two from copper- 
nickel, two from aluminum and one from stainless steel. 


Each tube had a 133.4 mm length exposed to steam. Tubes 45, 
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46, 47, 52 thru 56 were finned tubes with bore diameter of 
12.7 mm and fin root diameter of a 19.05 mm. Tubes 49, 51, 
and 57 thru 62 were finned tubes with fin root diameter of 
13.7 mm. Table II lists all the wire-wrapped tubes tested 
and their dimension. Tubes 63 thru 71 had smooth exteriors 
with a bore diameter of 12.7 mm and an outside diameter of a 
19.05 mm. Tube 36 and tubes 45 through 47 consisted of a 
family of spiral tubes, with triangular-shaped fins. 
Figure 3.8 shows a photograph of these four tubes which had 
a fin height of 1.0 mm and pitch of 1.06, 1.6, 2.1, and 2.5 
mm (tubes 36, 45, 46, and 47). 

Two commercially available tubes, manufactured by High 
Performance Tubes, Inc., were tested to investigate their 
heat-transfer performance compared to the other machined 
finned tubes. They were finned tubes (tubes 49 and 51) with 
fin height of a 0.75 mm and fin density of 1102 and 1181 
fins/m, respectively. А smooth tube (tube 50) was also 
Prepared by machining off the fins so that the effect of 
fins can be determined. These three tubes had the same 
outside root diameter (17.5 mm). In addition to these 
tubes, to study the effect of fin shape on the heat-transfer 
performance, four finned tubes were manufactured to complete 
two sets of tubes with different fin profiles. The first 
set consists of the tubes 06, 54, 55, and 56 with rectan- 
gular, parabolic, trapezoidal, and triangular fin-shapes, 
respectively, with a fin-base thickness of 1.0 mm. The 
other set consists of the tubes 17, 38, 53, and 52 with 
rectangular, parabolic, trapezoidal, and triangular  fin- 
shapes, respectively, with a fin-base thickness of 0.5 mm. 
Each tube in these two sets has the same fin height of 
1.0 mm and fin base spacing of 1.5 mm. Figures 3.4 shows a 
photograph of tubes 06, 54, 55, 56, while Figure 3.5 show 
schematic cross-sectional of the "parabolic" fin. Figures 


3.6 and 3.7 show a photograph and cross-sectional views of 
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tubes 17, 38, 52, and 53, respectively. As can be seen from 
Figure 3.5 and 3.7 these tubes do not have the exact shapes 
as stated above. For example, careful examination of the 
schematic cross-sectional of the parabolic tube (tube 54) 
showed that the fin-shape is almost straight near the fin 
base and circular at the fin tip, while the parabolic fin in 
Figure 3.7 shows that it had almost straight sides near the 
fin base with a sharp leading edge. Also, as seen in the 
same figure, a distinction between triangular and trape- 
zoidal fins is not possible. The reason for this is that 
the very thin fins lead to nearly the same fin thickness at 
their tips, because of the difficulties associated with the 
machining process. 

In order to test the effect of fin-metal thermal conduc- 
tivity on the heat-transfer performance, four spiral tubes 
with triangular fin profile were manufactured: one each from 
copper, copper-nickel , stainless steel and aluminum (tubes 
57, 58, 59 and 60, respectively). In addition, two tubes 
with rectangular fin profiles from  copper-nickel and 
aluminum were manufactured (Tubes 61 and 62, respectively). 
Due to the low thermal conductivity of these tubes, a 
smaller outside diameter (13.5 mm) was selected to minimize 
the tube metal resistance (note that these tubes have the 
Same nominal inside diameter as the other tubes). The 
Spiral tubes had a fin base thickness of 2.1 mm, while the 
rectangular fins had a fin thickness of 1.0 mm and a fin 
spacing of 1.5 mm. Figure 3.8 shows a photograph of these 
tubes. 

Finally, nine smooth copper tubes were manufactured, and 
each was wrapped with a titanium wire (0.5, 1.0 ог 1.6 mm 
diameter) at a nominal wire spacing of 1.0, 2.0 or 3.0 mm. 
Photographs of these wire-wrapped tubes are shown in Figures 
3.9, 3.10 Капа Еи 
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Photographs of Tubes with Different Fin 


Figure 3.4 


(tg - 1.0 mm). 


Shapes 


Figure 3.5 





Tracing of the Fin Profile of Tube 54. 
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Photographs of Tubes with 
Shapes (ty - 0.5 mm). 


Different Fin 


Figure 3.6 





(d) Trapezoidal (Tube 53) 


Figure 3.7 Cross-Sectional Photographs of Tubes with 


Different Fin Shapes Ср = 0.5 mm). 
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Photograph of Tubes with Spiral 


Triangular Fins (e - 1.0 mm). 


Figure 3.8 
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Photograph of Wire-Wrapped Tubes 


Figure 3.11 


(D. = 1.6 mm). 


IV. SYSTEM OPERATION AND DATA REDUCTION 


A. SYSTEM OPERATION 

Since copper has poor wetting characteristics with 
water, steam will normally condense on copper under a 
partial dropwise condensation mode, which is more effective 
than the filmwise condensation mode. Since the purpose of 
this investigation was to take data with filmwise condensa- 
tion and most of the tested tubes were made from copper, 
great care had to be taken to ensure that the filmwise 
condensation mode was in fact occurring. In order to ensure 
this, the tubes had to be treated according to the following 
procedure: 


1. The tube was rinsed with tap water to remove any 
contaminants that are soluble in the water. 


2. A mixture of equal parts of sodium-hydroxide апа есһу1 
alcohol was prepared and heated to about 80°C while 
frequently being stirred until it became watery. 


3. AÀ coating of this mixture was applied uniformly 
around the tube. 


4. The tube was placed in a steam bath and was heated by 
the steam for about an hour. 


5. А new coating was applied to the tube every 10 
minutes. 


6. The tube was then rinsed with distilled water and put 
immediately into the test section to avoid any contam- 
inants depositing on the tube which may lead to the 
dropwise problem. 

This process resulted in the formation of a thin layer of 
dark oxide that has high wetting characteristics. Since 
this layer was thin, its thermal resistance was negligible. 
This procedure was followed each time prior to the installa- 
tion of a tube. However, when re-testing ап already- 
darkened tube, is was heated in the steam only for 15 to 20 
minutes. 

Following the procedure described by Georgiadis [5], the 


test apparatus was brought to operating pressure апа 
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temperature by adjusting the input power to the boiler 
heaters, the cooling water flow rate through the tube, and 
the cooling water flow rate to the auxiliary condenser. 
Steady-state conditions were assumed when the operating 
conditions were stabilized with a steam temperature 
variation of + 2 yw V and a temperature rise ог the cooling 
water of + 0.005 K and + 0.01 K for atmospheric and vacuun, 
respectively. Once steady-state operating conditions were 
reached, the cooling water flow rate through the test tube 
was fed to the computer manually while the temperature rise 
of the cooling water through the test tube, vapor pressure 
and temperature were gathered automatically by the data 
aquisition system. For cooling water flow rates of 80 
percent (4.44 m/s for 19 mm O.D. and 12.7 mm I.D. tubes, and 
2.84 m/s for 17.5 mm O.D. with 15.6 mm I.D. tubes), 70, 62, 
54, 45,-35, 26, and 20 percent, and again 80 percent, two 
sets of data were taken. These cooling water flow rates 
were selected to give approximately equally-spaced heat flux 
values. After each change of the cooling water flow rate 
through the tube, the system pressure experienced a slow 
drift; so an adjustment of the water flow rate through the 
auxiliary condenser was required to maintain the system 
pressure at the operating pressure. AS mentioned in Chapter 
III, a view port was provided for visual observation to 
ensure filmwise condensation. Before each data collection, 
the appearance of the film was checked. If the film 
appeared to be patchy or there was an indication of dropwise 
condensation, the run was discontinued and the data were 
discarded. However, there were cases where the film 
appeared filmwise but the data collected at the end of the 
run (cooling water flow rate of 80 percent) were different 
from that collected at the beginning for the same flow rate. 
For example, the heat-transfer coefficient was as much as 10 


percent greater for the last data point and the cooling 
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water temperature rise was also greater than that measured 
at the beginning. As discussed by Georgiadis [5], this 
increased heat-transfer coefficient appears to be a result 
of the tube undergoing partial dropwise condensationwith 
exposure to the steam. Since this trend was observed even 
though no droples were visible, it is possible that the 
dropwise condensation was taking place at a microscopic 
level, especially near fin edges with a very thin condensate 
film. This phenomenon was observed primarily for the runs 
which followed the first treatment for darkening, probably 
due to the contaminants of the machine shop and since not a 
good layer of the dark oxide was obtained. All data 
presented in this thesis displayed less than 3% disagreement 
in the steam-side heat-transfer coefficient between initial 


айа final data sets. 


B. DATA REDUCTION | 

Initially, the program used for data reduction was the 
same as that used by Flook [6] including property functions, 
calibration curves for the cooling water flowmeter and for 
all thermocouples as well as the temperature rise due to 
frictional heating within the mixing chamber. Since tubes 
49, 50 and 51 had a different inside diameter than the tubes 
tested previously, and tubes 57 thru 62 were manufactured 
from metals with different thermal conductivity, the program 
was modified to include options for different tube diameter, 
thermal conductivity, fin shapes, and Sieder-Tate constant 
for the inside heat-transfer coefficient. 

The separation of the individual thermal resistances 
(water-side, wall, and vapor-side) from the overall heat- 
transfer resistance is very important in order to obtain 
expressions for the vapor-side heat-transfer coefficient. 
The overall heat-transfer resistance is given by equation 
(4.1), while the inside heat-transfer coefficient is given 


by a Sieder-Tate type equation (equation (4.2)). The value 


73 


of the leading coefficient C; must be known in order to 
calculate the inside heat-transfer coefficient from equation 
(4.2) and consequently the outside heat-tranfer coefficient. 
In order to determine the value of Ci, two methods were 
considered: the "direct" method and the "modified Wilson 
plot" method. 

















h; D 0. 15 
k y 


1. "Direct" Method 

This method is used to find the leading coefficient 
for the Sieder-Tate equation from an instrumented tube. As 
described by Georgiadis [5], a thick-wall smooth tube was 
manufactured with six thermocouples inserted into channels 
around the periphery of the tube. The average wall tempera- 
ture was found by averaging the temperatures indicated by 
the six thermocouples. He showed that the wall temperature 
distribution followed a cosine curve given by equation (4.3) 
with a maximum drop of 18 K between the top and bottom of 
the tube. 





AT қ со58 (4.3) 
Ат ит 
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where a is found to be from 0.135 to 0.202 and from 0.115 to 
0.179 under vacuum and at atmospheric pressure, respec- 
tively, in order to fit the temperature measurements. Using 
the average reading of all the thermocouples, Georgiadis [5] 
found values of 0.0635 and 26.4 for the leading Sieder-Tate 
coefficient and the constant B in equation (4.2), respec- 
tively, based on two runs each under vacuum and at atmos- 
pheric pressure. The value of C (0.0635) is greater than 
the well-known Sieder-Tate constant of 0.027 for the plain 


tubes, mainly owing to the coiled insert. The constant 
B = 26.4 is used for improved fitting of the experimental 
data. 


2. "Modified Wilson Plot" Method 
This method is a modification of the original Wilson 
plot method as modified by Briggs and Young [33] to accept 
data collected at various flow rates and temperatures. A 
Sieder-Tate equation was used for the inside heat-transfer 
coefficient, while a Nusselt type equation was used for the 


outside heat-transfer coefficient as given by equation 


(4.4): 


_ кр оғ (оғ - ру) Вер DE А 
ho = s ығ Dg q | | ) 


Both constants in equations (4.2) and (4.4) had to be deter- 
mined iteratively. Substituting equation (4.2) (with B 
- 0.0) and equation (4.4) in the equation for the overall 
heat transfer resistance given by equation (4.1), results in 


equation (4.5) below: 
аша НЫСЫН Lr (4.5) 
(T, кы) т Cy ke 2 И В 
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Equation (4.5) is a linear equation of the form: 





Ү = mX +b (4.8) 
where 
| . 
Y nc сет (4.9) 
бо Г 
“ &.10 
Я ( ) 
1 
т ЕС (4.11) 
y 


To begin the iteration, reasonable values меге 


assumed for Ci and f. With these values, the Y and X values 
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were calculated and a least-square technique was used to 
compute the slope and intercept values in equation (4.8). A 
new set of values for Ci; and f was then computed according 
to equations (4.11) and (4.12). This procedure was repeated 
until the assumed and computed values for Ci and f agreed 
to within 0.1 percent. Based on the Nusselt theory, f 
should take a value of 0.655, which is true for the case of 
zero vapor shear. However, for the experimental conditions 
in this thesis, the vapor velocity was from 1l m/s (atmos- 
pheric runs) to 2 m/s (vacuum runs), thus resulting in p 
values as high as 0.75. To account for the vapor shear 
properly, a correlation developed by Fujii and Honda [34] 
was also considered as shown below: 


_ Nu И Кашу ° — 
Retp е ` " ( ы ) 


š 


Using a Nusselt-type equation for the steam-side coeffi- 
cient, Georgiadis [5] found the leading coefficient (for 
equation (4.2)) С; to be 0.071, with the B value set equal 
to zero for a smooth tube. This C; value resulted in an 
inside heat transfer coefficient up to 6 percent greater 
than that based on the direct method.  Flook later used a 
Fujii-type [34] equation instead of the Nusselt-type equa- 
tion for the steam-side coefficient. This resulted in a 
slightly higher value (up to 3 percent) for the leading 
coefficient C;. The program used for the data collection 
allows an option for selecting either the Fujii-type or 
Nusselt-type equation for the steam-side coefficient. 
Georgiadis [5] and Flook [6] thought that the the 
"direct" method is more reliable, so the values of 0.0635 


and 26.4 were used for the constants С. апа В, 
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respectively, for the data reduction for tubes with an 
inside diameter of 12.7 mm. 

Later, Flook [6] used the "modified Wilson plot" 
method, with the Fujii-type equation for the steam-side 
coefficient to find the leading coefficient for the 
Sieder-tate equation for a copper tube and a stainless steel 
tube with thin (i.e., 0.5 mm thickness) tube walls since it 
was not possible to manufacture an instrumented tube for 
using the "direct" method to find the leading coefficient 
for the equation (4.1). For the copper and stainless steel 
tubes tube, values of C; = 0.0756 and СЕ 0.0688 меге 
obtained respectively. For this thesis effort, initially 
the values of C; = 0.0635 and B = 24.6 were found with the 
"direct method" as said above, but finally the "Modified 
Wilson Plot" directly on the finned tubes was used. This 
was the same method as described above but for each tube a 
different Sieder-Tate coefficient was found and used for the 


data reduction. 
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V. RESULTS AND DISCUSSION 


А. INTRODUCTION 


During this thesis effort, a number of data runs were 
made using the procedures described in Chapter IV. Each 
tube was tested at least three times, both under vacuum and 
at atmospheric pressure on different days, to ensure repeat- 
ability of the data. Complete filmwise condensation condi- 
tions were maintained, and the non-condensing gas 
concentration was calculated at the beginning and at the end 
of each run for every tube to ensure there were no major 
leaks in the system. The computed mass concentration of 
non-condensing gases was kept between O0 and - 2.5 percent in 
order for the data run to be accepted. The mass concentra- 
tion of the non-condensing gases was always a negative 
number because of the slight inaccuracies in the measurement 
of pressure and steam temperature. As discussed earlier in 
Chapter III, the test apparatus would allow only a negli- 
gible amount of non-condensing gases to be leaked into the 
apparatus. Since continuous venting was provided throughout 
all runs (see Chapter III), build up of non-condensing gases 
was not possible. 

A summary of finned tubes tested by Georgiadis [5] and 
Flook [6] and those tested during this thesis effort, as 
well as the resulting enhancements are provided in Table 
TIT. Further, Table IV presents the wire-wrapped tubes 


tested and their heat-transfer performance. 
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B. WATER-SIDE HEAT-TRANSFER COEFFICIENTS 

As mentioned in Chapter IV, Georgiadis [5] used two 
methods to find the Sieder-Tate coefficient: the "direct" 
method and the "modified Wilson plot" method The "direct" 
method involved the measurement of the average tube wall 
temperature using six thermocouples embedded within the wall 
of a smooth tube. He took data with filmwise condensation 
occurring outside. He changed the water velocity from about 
0.8 to 4.5 m/s and correlated the data resulting in a 
Sieder-Tate constant C; of 0.064 (see equation (4.2)) with a 
B value of 26.4. Also, taking data on an unistrumented 
smooth tube, and making a modified Wilson analysis, he found 
a Sieder-Tate constant of 0.071 with the B value set equal 
to zero. When the h; values computed using the results of 
these two methods were compared, they agreed to within 6 
percent. 

During the present study, however, a.third approach was 
tested. For this purpose, the modified Wilson analysis was 
carried out directly on finned-tube data, resulting in E 
values around 0.069 with the B value set equal to zero. The 
hi values computed using this analysis in fact lie between 
the values computed by the two methods decribed earlier. 
For this reason, the third method was used throughout this 
investigation іп computing the outside  heat-transfer 
coefficient. 

Since the water-side geometry for tubes 45 through 47, 
49 through 51, 52 through 56 , 57 through 62 and 63 through 
71 was essentially the same, the values of Ci should be very 
nearly the same for all these tubes. However, circumferen- 
tial temperature variations can influence the inside coeffi- 
cient, and these variations will depend on the condensate 
retention angle. As discussed in Chapter II, the retention 
angle is strictly dependent on the fin spacing. Therefore, 


the experimentally found С; values should vary from tube to 
tube. 
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C. REPEATABILITY OF DTS 

In order to insure the reliability of the data taken, 
all data runs were repeated, as Georgiadis [5] and Flook 
[6] did, at least three times on different days. The 
computed steam-side coefficients for similar conditions 
(i.e., same tube and about the same operating conditions) on 
different days agreed to within + 5 percent for some tubes 
and + 10 percent for others. Additionally, data runs were 
performed on six finned tubes (tubes 6, 17, 27, 28, 36, 38) 
under similar conditions to verify the repeatability with 
data taken by Flook [6]. Georgiadis [5] also had tested 
tubes 6 and 17. Georgiadis and Flook processed their data 
using the Sieder-Tate constant found by the "direct" method 
(see Section B of Chapter IV). Therefore, in order to 
perform.a fair comparison, their data were reprocessed by 
the method used during this investigation (i.e., using the 
"modified Wilson plot" method directly on finned tube data). 
Figure 5.1 shows the experimental steam-side heat-transfer 
coefficients of Georgiadis [5] and Flook [6] and those 
obtained during this investigation for tube number 6 under 
vacuum. Figure 5.2 shows similar data for tube 6 under 
atmospheric pressure, whereas Figures 5.3 and 5.4 show 
Similar data for tube 17. For comparison purposes, the 
smooth tube data were plotted as well as the data predicted 
by Nusselt theory. As can be seen from Figure 5.1, the data 
obtained during this work under vacuum conditions fall above 
and below the data of Georgiadis [5] and Flook [6] with a 
maximum variation of 20 percent. At atmospheric conditions, 
the data fall about 15-20 percent below those of Georgiadis 
[5]. The agreement with tube 17 is much better as can be 
seen from Figures 5.3 and 5.4. For tubes 27, 28, 36, and 38 
(see Figures 5.5 through 5.9) the agreement is very good, 


with a deviation of only + 5 percent except for tube number 
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ГА 


36 at atmospheric pressure. In this case, the agreement was 
+ 26 percent. The disagreement, especially for tube 6, is 
probably due to a very small contamination on the outside 
tube surface area of the tube leading to partial dropwise 


condensation conditions during the runs made by Georgiadis 


and Flook. Typical uncertainty bands are also included at 
low and high heat flux to indicate the maximum possible 
uncertainty. However, the repeatability for most of the 


tubes shows that the experimental uncertainty is always less 
than that indicated by the uncertainty bands. 

In these figures, the steam-side coefficient is plotted 
versus the heat flux, and as is always the case for conden- 
sation, the heat-transfer coefficient decreases as the heat 
flux increases. The curves shown in these figures (and 
subsequent figures) are the  least-squares-fit curves 
according to the following equation: 


- 


а = а ATP (S) 


š 


where AT was computed using the following equation: 


q = ho AT (5227) 


where q was measured experimentally and the steam-side coef- 
ficient was calculated by subtracting the inside and wall 
thermal resistances from the overall resistance as given by 


equation (5.3): 


1 
һ = z o 5 ad 
о Leg. Sad ( ) 
where Uo ТЕ 
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Ry = e ori (5.4) 
In equation (5.1) the coefficient a and the exponent b are 
experimentally determined constants. The values of a and b 
both under vacuum and atmospheric conditions for all the 
finned tubes and for the wire-wrapped tubes tested are given 
in Tables V and VI, respectively. 


TABLE V 
CONSTANTS OF EQUATION (5.1) FOR FINNED TUBES TESTED 





Tube VACUUM RUNS ATMOSPHERIC RUNS 

Tube S t e 

No. (mm) (mm) (mm) 
06 | 1.50 1.00 1.00 61190 0.72588 91685 0: 7200 
17 | 1.50 0.50 1.00 63210 0.72%65 95079 0.71433 
37 | 1.50 1.00 0.50 49315 DN I2 wu TUE 
28: 52595 1.00 0.50 52515 o77 wearer 
36 | 0.00 2.10 1.00 63847 0.72346 93246 0.71243 
38 | 1.50 0.50 1.00 67881 0.72661 | 102870 0.71715 
4S | 0.00 2.50 1.00 63459 0.22502 99843 0.70538 
46 | 0.00 1.60 1.00 65976 0.22740 | 101200 0.71728 
ІІ | 70200 16 1.00 63456 0.72683 98200 0.71167 
49 | 0.51 0.34 1.00 61325 0.72920 89327 0.71046 
50 | ----- | -- --- | ----- 28144 0.70192 43092 0.661 
81 | 0.59 0.32 1.00 61928 0.72950 93058 0.71267 
S” | 50 0.50 1.00 63672 0.72508 95028 0.71388 
s3 | 1 50 0.50 1.00 62772 0.72613 94397 O 71270 
54 | 1.50 1.00 1.00 61452 0.72125 94299 0.71061 
55.11.50 1.00 1.00 66987 0.72326 90906 0.70660 
56 | 1.50 1.00 1.00 65674 0.72303 94829 0.70866 
57 | 0.00 2:10 1.00 65871 0.72540 94095 0.69296 
58 | 0.00 2 10 1.00 52916 0.70566 nass 0.68875 
59 | 0.00 2.10 1.00 25364 0767315 38902 0.67594 
60 | 1.50 1.00 1.00 57695 0.71805 82971 0.69082 
61 | 1.50 1.00 1.00 40173 0.70228 72273 0. 69034 
62 | 1.50 1.00 1.00 50021 0.71093 85279 0.69645 


—— 00 mn mm eae SE h 


D. EFFECT OF FIN PITCH ON HEAT TRANSFER PERFORMANCE OF 
SPIRAL TUBES WITH TRIANGULAR-SHAPED FINS 


This section presents results showing the variations of 
the steam-side heat-transfer coefficient with heat flux 


having fin pitch as a parameter. Data were taken on four 
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TABLE VI 


CONSTANTS OF EM 


FOR WIRE-WRAPPED TUBES P 


ESTED 


VACUUM RUNS ATMOSPHERIC RUNS 
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63 1-6 0. 2.54 26207 0.68301 86280 0.65983 
64 16 2 3.63 32515 0.69452 322522 0.67158 
65 1-6 3. 8.62 32815 0.69733 52403 0.66672 
66 1.0 0. 1.25 27097 0.68988 47762 065921 
67 1.0 1. 2.82 35403 0.70265 58145 0.67488 
68 1.0 2, 3.91 20521 0.70144 58421 0.67219 
69 0.5 r 1.60 34291 0.59225 59000 0.66327 
70 0.5 B. 2.54 40254 0.70477 63569 0.67866 
p 0.5 a 3.63 39179 0.70190 62272 0.67547 








copper tubes with spiral triangular fins. These tubes have 


the same fin height of 1.0 mm and have fin pitches of 1.06, 
1.6, 2.1, and 2.5 mn. Figures 5.10 and 5.11 present data 


for these tubes under vacuum and at atmospheric pressure, 


respectively. The smooth-tube data and a curve representing 


Nusselt theory also included for comparison. The best heat- 
transfer performance was obtained with the tube with a fin 
pitch of 1.6 mm. 


As shown in Appendix C, the uncertainty in the calcula- 


tion of the steam-side coefficient increases as the heat 


flux decreases. Therefore, the comparison of the finned 


tubes should be performed at a high heat flux, where the 
uncertainty is small. 
through the 


the 


The comparison of finned tubes is 
This 


steam-Side heat-transfer coefficient of a 


made enhancement ratio, Eo. ratio is 


defined as 
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finned tube to that of the smooth tube (same diameter as the 
finned tube root diameter) at the same heat flux. Heat flux 
values of 0.25 and 0.75 MW / m2 were chosen for vacuum and at 
atmospheric conditions, respectively. For the spiral trian- 
gular fins, maximum enhancement ratios of about 3.9 and 6.1 
under vacuum and at atmospheric pressure were found. The 
enhancement ratio at atmosheric pressure is always higher 
than that under vacuum conditions. At atmospheric pressure, 
a higher temperature exists, so the condensate has a smaller 
viscosity, which results in improved drainage from the fin 
valleys, and smaller surface tension which results in a 
smaller retention angle. As discussed in Chapter II, the 
flooded portion of the tube has another thermal resistance 
due to the thick layer of condensate. Reducing the flooded 
portion of the tube increases the heat-transfer performance. 

Cross plots of the enhancement ratio versus fin pitch 
are shown in Figure 5.12, while Figure 5.13 shows a cross 
plot of the normalized ratio Eo/Ar (the ratio of the 
enhancement ratio to the area ratio). Generally, as the fin 
pitch increases, Eo/Ar increases. Table III shows that tube 
47, with a fin pitch of 1.06 mm, has the largest area ratio, 
while tube 45 with fin pitch of 2.5 mm has the smallest area 
ratio, and tube 47 has a poorer performance than either 
tubes 45 or 46. The poor performance shown by tube 47 can 
be explained by the fact that, as the pitch decreases, the 
area of the tube increases, but at the same time the reten- 
tion angle increases and more flooding occurs. This means 
that as the fin pitch increases, the interfin spacing was 
covered by a thick layer of condensate. The additional 
thermal resistance induced by this layer of condensate over- 
powers the benefit gained from increased surface area, so 
the heat-transfer performance is reduced. 

As the pitch increases from 1.05 mm to 1.6 mm, the 


retention angle decreases more than the area decreases. 
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This results in a larger enhancement ratio for the tube with 
a pitch of 1.6 mm. Beyond this point, the area ratio 
decreases while the retention angle decreases slowly and 
results in a smaller heat transfer performance. In order to 
obtain a clearer understanding on the heat-transfer perform- 
ance, the enhancement ratio was divided by the area ratio, 
and thus, the effect of the changing area was eliminated. 
As Figure 5.13 shows, the enhancement ratio depends on other 
factors in addition to the fin area, such as the surface 
tension effect. As discussed in Chapter IL, the surface 
tension-induced pressure gradient from the fin tip to the 
fin root is responsible for thinning of the condensate film 
and thereby improving the heat-transfer performance in the 
unflooded portion of the tube. Also, the surface-tension 
forces cause condensate flooding, resulting in poorer 
performance in the flooded portion of the tube. Figure 5.13 
shows that the normalized enhancement ratio is higher for 
the tube with fin pitches of 2.1 mm and 2.5 mm under vacuum 
and at atmospheric pressure, respectively. Therefore, the 
optimum fin pitch is between 2.1 and 2.5 mm based on normal- 
ized enhancement ratio, while the optimum fin pitch is 
1.6 mm based on the enhancement ratio. As shown by Edwards 
et al. [20], as the pitch increases for the same fin height, 
the heat-transfer coefficient increases. However, as the 
fin pitch increases, the tube is easily flooded. Because the 
retention angle is greater under vacuum than that at atmos- 
pheric pressure, the tube with fin pitch of 2.5 mm has a 
smaller retention angle than that of the tube with fin pitch 
of 2.1 mm under vacuum, resulting in poorer heat-transfer 
performance. However, at atmospheric pressure, since the 
retention angle is less than under vacuum conditions, the 
tube with a fin pitch of 2.5 mm has a better heat-transfer 


performance than the tube with a fin pitch of 2.1 mm. 
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E. EFFECT OF FIN SHAPE ON HEAT-TRANSFER PERFORMANCE 

In order to study the effect of fin shape, data were 
taken on two sets of copper tubes with fins of four 
different shapes. All the fins were manufactured with same 
fin spacing and fin height. The first set of tubes (17, 38, 
52 and 53) had rectangular, parabolic, triangular and trape- 
zoidal fin shapes, respectively, with a fin-base thickness 
of 0.5 mm, while the second set of tubes (6, 54, 55 and 56) 
had a fin-base thickness of 1.0 mm. Dimensions for these 
fins are given in Table III. ‚ 

The performance of tubes 17, 38, 52, and 53 under vacuum 
conditions is shown in Figure 5.14, while Figure 5.15 
depicts their performance at atmospheric pressure. For 
comparison purposes, data for a smooth tube are also shown. 
The tube with the "parabolic" fin profile (tube 38) showed 
the best heat.transfer performance. while the other three 
tubes performed about equally, under both pressure condi- 
tions. As shown in Table III, an enhancement ratio of 4.1 
and 6.2 were obtained for the tube with "parabolic" fins 
under vacuum and at atmospheric pressure, respectively. 
Also, it can be seen that the area ratio of the tube with 
"Parabolic" fins is less than that of the tube with rectan- 
gular fins and the tube with trapezoidal fins, and larger 
than that of the tube with triangular fins. Eliminating the 
effect of increased area, the ratio Eo/Ar is larger for the 
tube with "parabolic" fins under atmospheric and vacuum 
conditions. 

The reason for the greater enhancenent ratio is probably 
due to the continuous change of radius of curvature 
(increasing from the fin tip to the fin root) for the "рака- 
bolic" shaped fins. The condensate film has a convex shape 
at the fin tip and a concave shape at the fin root. The 
condensate film follows approximately the curvature of the 
wall surface at the fin tip. Because of the convex conden- 


sate surface at the fin tip and the concave condensate 
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surface near the fin root, the pressure within the conden- 
sate is larger at the fin tip and smaller at the fin root. 
Therefore, an appreciable pressure gradient exists from the 
fin tip to the fin root. The gradual increase of radius of 
curvature results in a gradual decrease in pressure within 
the condensate, which is very important for improved conden- 
sate flow, resulting in a thinner film and larger heat- 
transfer coefficient than if the fin sides were flat. 
Therefore, the parabolic fins should outperform all other 
three tubes, as shown in Figures 5.14 and 5.15. 

Also, Adamek [24] and Mori et al. [8] have shown that 
the optimum fin shape is that which induces a continuous 
pressure gradient due to the surface tension effect or which 
has large curvature at the fin tip and continuously 
decreasing toward the fin root. Therefore a continous 
decrease in the pressure gradient exists and this thins the 
condensate film continuously, resulting in better heat 
transfer performance. For the case of tube 38, there is a 
continuous decrease of the curvature, while this was not 
happening for the other three tubes of the first set 
although they had a sharp leading edge. Therefore tube 38 
exhibits better heat-transfer performance than the other 
three tubes in its group. 

Figures 5.16 and 5.17 show the performance of tubes 6, 
54, 55 and 56 under vacuum and at atmospheric conditions, 
respectively. Again, the smooth tube data are included for 
comparison purposes. For this set of tubes, the best 
performance is obtained from tube 55 with the trapezoidal 
fins, while the tube with triangular fins performed second 
and the remaining two tubes performed about equally well 
under vacuum conditions. However, the parabolic and trian- 
gular shapes outperformed the other shapes for atmospheric 
conditions. The poorer performance of tube 54 under vacuum 


conditions was not expected. As Figure 3.7 shows, the fins 
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do not have the shapes as claimed above due to the difficul- 
ties encountered in machining. The unexplainable trends 
shown by the tube with "parabolic" fins (tube 54) was found 
to be the result of the actual fin shape that was very 
different from what was expected. At the conlusions of the 
data runs presented in this thesis, this tube was destroyed 
and a magnified photogragh of the fin cross-section was 
taken. AS can be seen from Figure 3.5, these fins do not 
have a profile with gradually decreasing curvature from the 
fin tip to root. Therefore, unlike in the previous set of 
tubes (17, 38, 52 and 53), the data taken on the second set 
of tubes are inconclusive. Since the condensation process 
on a finned tube is extremely complex, owing to the very 
large number of parameters, it may be unwise to draw conclu- 
sions from the above-mentioned results. Nevertheless, the 
data for the second set of tubes are presented in this 
thesis for completeness. | 

From Figure 3.7 it is clear that: 1) the "parabolic" 
fins of the tube 38 with fin base thickness of 0.5 mm had a 
straight fin side and a fin tip with a sharp leading edge, 
while careful examination of a cross-section of the "рака- 
bolic" fins of tube 54 shows a straight fin side with a 
near semicircle at the fin tip. Tube 38 had a small radius 
of curvature at the fin tip, while tube 54 had larger 
radius of curvature at the fin tip. These differences in 


geometry may have caused the observed data. 


F. EFFECT OF FIN THERMAL CONDUCTIVITY ON PERFORMANCE 

To investigate the effect of fin-metal thermal conduc- 
tivity on the heat-transfer performance, four spirally 
finned tubes with triangular fins and two tubes with rectan- 
gular fins were manufactured. As shown in Table I, the 
four spiral tubes were made of copper,  copper-nickel, 
Stainless steel, and aluminum, respectively (tubes 57, 58, 
59 and 60). The two tubes with rectangular fins were made 


of copper-nickel and aluminum (tubes 61 and 62). All tubes 
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had the same fin height of 1.0 mm. Тһе tubes with a rectan- 
gular fin shape have the same fin dimensions as the 
"optimum" copper tube (tube 6) found by Georgiadis (51, 
while the spiral tubes had a fin pitch of 2.1 mm. The 
results for data runs taken under vacuum and at atmospheric 
pressure are shown in Figures 5.18 and 5.19, respectively, 
for the spiral tubes. For the tubes with rectangular fins, 
the variation of heat-transfer coefficient with heat flux is 
shown in Figures 5.20 and 5.21 under vacuum and at atmos- 
pheric conditions, respectively. Figure 5.18 shows that the 
copper tube exhibits the best heat-transfer performance, 
while the stainless steel shows the worst performance. The 
second best is the aluminum spiral tube followed by the 
copper-nickel tube. The same trend is also seen in Figure 
5.19 at atmospheric pressure. Enhancement ratios as high 
as 3.5 and 4.4 under vacuum and atmospheric conditions, 
respectively, were found. Cross plots of enhancement ratio 
Eo and normalized enhancement ratio Eo/Ar versus the thermal 
conductivity аге shown іп Figures 5522 and 512 
respectively. 

A similar trend exists for tubes with rectangular fin 
profiles as shown in Figures 5.20 and 5.21. For comparison, 
the data of Flook [6] for the copper tube with D, - 14.5 mm 
and D; = 13.5 mm with the same fin dimensions as the 
aluminum and copper-nickel tube (tube 39) are also included. 
Also in the same Figures the data of tube 6 are included. 

The thermal conductivity of aluminum (167 W/m.K) is 
about half of copper (385 W/m.K), while copper nickel and 


stainless steel have much lower values (i.e., 45 W/m.K and 
15 W/m.K, respectively). Since the thermal resistance 
through the fin increases  (i.e., the fin efficiency 


decreases) with decreasing thermal conductivity, the copper 
tube must show the best heat-transfer performance, while the 


stainless steel tube must show the poorest performance. As 
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can be seen from Figures 5.18 through 5.22, this trend is 


very clear. 


G. PERFORMANCE OF WIRE WRAPPED TUBES 
1. Condensate Retention Angle for Wire-Wrapped Tubes 

Wire-wrapped tubes are somewhat similar to finned 
tubes with regard to their susceptibility to condensate 
retention. Since the portion of the tube with the retained 
condensate would result in a poor heat-transfer performance, 
1t was mecessary to study this phenomenon on these wire- 
wrapped tubes. For this purpose, it is possible to use the 
Webb et al. [17] model (as discussed in Chapter II), which 
was developed for a tube with fins of arbitrary shape (see 
equation (2.10)). However, they did not test this model for 
wire-wrapped tubes. Therefore, it was necessary to experi- 
mentally measure the condensate retention angle for these 
tubes, so that the Webb et al. model can be modified to 
predict it for wire-wrapped tubes. For this purpose, meas- 
urement of the condensate retention angles were made for the 
wire-wrapped tubes under static conditions. Photographs 
were taken; slides were made and accurate measurements were 
made on the screen. Figure 5.24 shows photographs of 
portions of tubes with a wire diameter of 0.5 mm and pitches 
of 1.6, 2.5 and 3.6 mm (tubes 69, 70 and 71), respectively. 
Careful examination of the photographs revealed that water 
was trapped all around the tube in a region very near of the 
wires. The film thickness was about the same along the 
distance from the top of the tube to the bottom, until some 
point where the meniscus of the water had a concave profile 
in the region where the tube starts to be fully flooded 
between two adjacent wires. In this case, the retention 
angle is defined as the angle from the bottom of the tube to 
point A (see Figure 5.24), where the surface tension forces 
balance the weight of the condensate. The retention angle 
was computed by equation (5.5) based on the height from the 
bottom of the tube to point A and the tube radius. 
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TABLE VII 
MEASURED RETENTION ANGLES (WV) 


Tube # Dw Р Measured Calculated 
# (mm) (mm) Degrees Degrees 

63 1.6 Das 116 

64 1.6 Jj 81 1% 

65 lk. 6 4.6 Z3 57 

66 1.0 2.0 118 123 

67 i. O0 228 87 78 

68 1.0 3059 70 29 

о HB n м do 

71 бе 5 3.6 47 54 

x 
کڪ‎ е z)? (5.5) 


Table VII lists the results for all nine tubes (tubes 63 
through 70). Using the Webb et al. model [14], attempts 
were made to predict the measured retention angles, and 
comparison is shown in Figure 5.25 (see the "triangular" 
symbols--the "star" symbols will be discussed below). As 
can be seen, the Webb et al. model overpredicts the conden- 
sate retention angle for most of the tubes. Therefore, this 
model was modified for the present study (i.e., steam 
condensation on horizontal wire-wrapped tubes). The exami- 
nation of equation (2.10), reveals that the condensate 
retention angle would decrease with increasing wire diameter 
while all other parameters are kept constant. However, the 
experimentally measured retention angles shown in Table VII 
show the opposite trend (i.e., y increases with increasing 
Dy). Therefore, it was necessary to modify the Webb et al. 


model as shown in equation (5.6). 


- S ONU МЕС) Dw) 2 5.6 
pre fine om (BY S| 9 
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where а апа [f are empirical constants to fit the measured 
retention angle. In order to compute these a and b values, 
a numerical procedure was followed, by minimizing the sum of 
squares of the deviations of the computed and measured y 
values. This procedure resulted in a and § values of 0.85 
and 0.18, respectively. The comparison between the values 
computed by equation (5.6) and the experimental values is 
shown in Figure 5.25. As can be seen, this equation agrees 
to within * 15 percent with the experimental data. 
2. Experimental Data of Wire-Wrapped Tubes 

This section presents results showing the variation 
of the steam-side coefficient with heat flux for the wire- 
wrapped tubes tested. Three wire diameters (0.5, 1.0 and 
1.6 mm) were used with three different pitches for each wire 
diameter. These are shown in Table IV; The tubes with a 
1.6 mm wire diameter had pitches of 2.5, 3.6, and 4.6 mm; 
the tubes with 1.0 mm wire diameter had pitches of 2.0, 2.8, 
and 3.9 mm; and, the tubes with 0.5 mm wire diameter had 
pitches of 1.6, 2.5, and 3.6 mm. Figures 5.26 and 5.27 show 
the variation of heat-transfer coefficient with heat flux 
under vacuum and at atmospheric pressure, respectively, for 
the tubes wrapped with a l.6-mm-diameter wire, while Figures 
5.28 and 5.29 show the variation of the heat-transfer coef- 
ficient under vacuum and at atmospheric conditions, respec- 
tively, for the tubes wrapped with 1.0-mm-wire diameter. 
Figures 5.30 and 5.31, show similar results for the tubes 
wrapped with 0.5-mm-diameter wire under vacuum and at atmos- 
pheric pressure, respectively. Figures 5.26 to 5.31 show 
that the best performance was obtained for the tube with a 
wire diameter of 0.5 mm and a pitch of 2.5 mm, both under 
vacuum and at atmospheric pressure. The second best perform- 


ance was given by the tube with a wire diameter of 0.5 mm 


EZO 


and a pitch of 3.6 mm under both pressure conditions. The 
third best performance was given by the tube with a wire 
diameter of 1.0 mm and a pitch of 2.8 mm under atmospheric 
pressure, and the tube with a wire diameter of 1.0 mm and a 
pitch of 3.9 mm, under vacuum pressure. The tubes with wire 
diameters of 1.6 mn and 1.0 mm with pitches of 2.5 and 
2.0 mm, respectively, showed the worst performance among all 
the tubes tested. 

As Figures 5.26 and 5.27 show, the performance of 
the tubes with a wire diameter of 1.6 mm and pitches of 3.6 


and 4.6 mm interchanged as the pressure conditions were 


changed from vacuum to near-atmospheric. The same trend 
happened for the tube with wire diameter of 1.0mm and 
pitches of 2.8, and 3.9 mm. This behavior can only be 


explained by the retention angle phenomenon. The tubes were 
wrapped in order to take advantage of the condensate thin- 
ning as a result of the varying surface-tension forces in 
the space between the wires. However, as mentioned in 
Chapter II, these surface-tension forces lead to a deleter- 
ious effect owing to condensate retention, especially when 
the fin spacing is small. If the extent of condensate 
retention is large, the enhancement gained over the 
unflooded portion of the tube (owing to condensate thinning) 
may be offset in poor performance in the flooded portion 
with retained condensate. As discussed in subsection 1 
above, the retention angle is higher under vacuum conditions 
for tubes with a pitch of 3.6 than for tubes with a pitch of 
4.6. This is true for tubes with pitches of 2.8 and 3.9 mm 
as well. However, under atmospheric pressure, the retention 
angle decreases resulting in less flooding, and better 
thermal performance. Figures 5.26 and 5.28 show that the 
tubes with wire diameters of 1.6 and 1.0 mm with pitches of 
2.5, and 2.0 mm, respectively, had worse performance than 


the smooth tube under vacuum. Again these are the results 
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of the surface tension effects. For these tubes the reten- 
tion angle was about 100 degrees, resulting in about half 
of the tube being effective for the heat transfer. The 
other half contributes a small amount to the heat transfer 
performance due to the thick layer of the condensate. The 
heat-transfer enhancement ratio is shown clearly in Figure 
5.32, as a function of p/d,. The optimum p/D,, appears to be 
near 5 to 6. 
3. Modifications to Fujii et al. Model 

As mentioned in Chapter II, Fujii et al. [32] devel- 
oped a semi-theoretical expression to predict succesfully 
their condensation data on wire-wrapped tubes, using ethanol 
and R-11 as the working fluids. Ав stated in Chapter II 
(Section E), the vapor-side enhancement could be easily 
computed using equation (2.65) provided A » 15. However, 
for most of the tubes tested during this investigation, the 
values of A were as low as 6. Thus, it was not possible to 
use equation (2.62) and equation (2.59) had to be used. 
When equation (2.59) is plotted as a function of A and фр, 
the result is shown in Figure (5.33). Notice that for 
А > 18, ШЕТ.) is independent of ф Also for A > 18 the 
functional dependence of F, on the angle ф сап be approxi- 


mated by a straight line: 


F1(6,A) * F1(0,4) - m à (5.7) 


where m is the slope of the straight line. In equation 
(5.7), F1(0,A) for all values of A is given by equation 
(226298 

Because of the flooding that occurs on the lower 


portion of the wire-wrapped tubes when steam is being 
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condensed, in order to compute the heat-transfer performance 
of the unflooded portion of the tube, the function Fj(09,A) 
must be integrated from ф = 0.0 Со ф = т - y, as needed for 
equation (5.8): 


l 
FACHA) = al” Fig) do (5.8) 


Notice that this equation is the same as equation (2.58) 
except that the integration is performed up to the angle @ 
(піса ist Br - we Fujii et al. neglected the retention- 
angle effect for their low-surface-tension fluids. As 
mentioned earlier, for A > 15, equation(5.7) is valid, but 
the slope m can depend on the value of A. In order to find 
the dependence, the slope was approximated for values of 
A > 5 and these slopes are plotted in Figure (5.34) as a 
function of A. Using a least square fit of the calculated 


data, the following functional form was derived: 


b | (323) 


where а = 0.177, апа 
b = -0.756 

Notice that the actual computation of m was performed only 
for 5 < A < 19. The numerical integration of F)(o,A) for A 
> 18 was not possible owing to overflow limitations of the 
computer. Further, the integration of Е] (ФА) Ғог А < 5 
was not performed since this was outside of the experimental 
conditions; the computed minimum A value was about 6. 
Figure 5.34 shows that the least-squares-fit curve has been 
extrapolated for A > 18. Even though such overextrapola- 


tions are not generally recommended, it appears reasonable 
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for this situation. The justification for this extapolation 
lies with the diminishing value of the slope m. 


Substituting equation (5.7) into (5.8) yields: 


к-2 


| 92 (SEPON) 
Е2(%А) х -1|Ғ1(0,4)9% - пз 1 


Ө 


In order to modify the Fujii et al. model, three assumptions 
were made: 1) heat transfer through the wire is negligible 
compared to that through the interwire space for the 
unflooded portion of the tube; 2) heat transfer through the 
flooded portion is by one-dimensionalconduction, and 3) in 
the flooded portion of the tube, heat transfer through the 
condensate between the wires is negligible compared to the 
heat flow across the wires. Based on assumption (1) above, 
the heat-transfer coefficient for the unflooded portion can 


be expressed as follows: 


fu S F2(0,A) (5 


— == =—ЭШЬ 


hg — (s + Dy) F2(0,0) 


where F4(0,0) is obtained numerically using equation (5.8). 
and where h. is the heat-transfer coefficient for the smooth 
tube predicted by Nusselt theory [27]. Also, the heat- 
transfer coefficient across the wires in the flooded portion 
of the tube can be computed by equation (5.12) (derived in 
Appendix A): 





1 Dy dx )5.12( 
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Now, combining equations (5.11) and equation (5.12),  equa- 
tion (5.13) can be obtained in order to express the average 


heat-transfer coefficient for the entire wire-wrapped tube. 





= = V 
h кел, (1 = p gle) = FD (5.13) 


Then the enhancement ratio is given by equation (5.14): 


hu V ће ф D, 





In order to fit the experimental data, a value of 0.02 was 
selected for the constant C in equation (2.61). Figure 5.35 
shows the experimental data and the values calculated from 
the modified Fujji et al. model. As can be seen from this 
figure, good agreement of the experimental data and the 
predicted values exists. However, a clear trend does not 
exist for the effect of wire diameter on the heat-transfer 
performance. The assumptions made during this study, and 
other possible mechanisms not taken into consideration owing 
to the complex nature of this problem, such as convective 


effects may be responsible for the observed comparison. 


H. HEAT-TRANSFER COEFFICIENT FOR COMMERCIAL TUBES 

As mentioned in Chapter IV, two commercially available 
finned tubes were tested. These tubes were manufactured by 
High Performance Tube, Inc., and had fin densities of 1.1 
and 1.18 fins/mm (tubes 51 and 49). Figures 5.36 and 5.37 
show the variation of the heat-transfer coefficient for 
these tubes with heat flux under vacuum and at atmospheric 
pressure, respectively. Data for a smooth tube (tube 50), 
are also shown. The best performance was given by the tube 


with fin density of 1.1 fins/mm under both pressure 
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conditions. Visual examination of the condensation process 
revealed that both finned tubes were fully flooded. This 
was also observed, under static conditions. Since these 
tubes had a high density of fins with small fin height, the 
surface tension effect is not important for the flow of the 
condensate along the fin side. However, the surface tension 
retained the condensate between the fins, so these tubes 
were fully flooded. Therefore, as Table III shows, although 
they have the largest area ratio, the heat-transfer 
performance is worse than most of the other copper finned 
tubes tested during this thesis effort. Nevertheless, even 
though these tubes were fully flooded, the normalized 
enhancement ratio (Eo/Ar) was greater than unity in both 
cases (especially for atmospheric pressure), indicating an 


enhancement greater than first due to an area increaSe. 
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А. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


The use of fins lead to significant enhancement of the 
steam-side heat-transfer coefficient, which is greater 
than, the area ratio (finned tube area to smooth tube 
area) despite condensate retention between fins. 


Enhancement ratios as high as 1.84 and 2.6 under 
vacuum and at atmospheric conditions, respectively, 
were realized for the wire-wrapped tube with a wire 
diameter of 0.5 mm and a ze of 5.1 тш. This rs due 
to the surface-tension.effect resulting in thinning of 
the condensate between wires. 


The Webb et al. [31] model was succesfully modified to 

predict the condensate retention angle for the 

aban rapped tubes. Maximum error of 15 percent was 
ound. 


The Fujii et al. [32] model was modified to predict 
the vapor-side  §heat-transfer coefficient for 
high-surface-tension fluids (i.e., water). A favorable 
eee berwecn the modified Fujii et al. model an 
the experimental data was found. 


The tube with a “parabolic” fin profile outperformed 
the tubes with күлала: trapezoidal and rectangular 
fin shapes. For the, tube with arabolic ins, 
enhancement ratios of 4.1 and 6.2 under vacuum and at 
atmospheric pressure, respectively, were obtained. 


For tubes with spiral triangular fins, the optimum 
pitch was found to be about of 1.6 mm. Enhancement 
ratios of 3.9 and, 6.1 under vacuum and at atmospheric 
pressure, respectively, were obtained for this tube. 


The High Performance tubes, although they have the 
highest area ratio among all the tubes tested, show 
poorer performance than most of the other tubes. This 
15 Шаа owing to the high condensate _ retention 
angle. nhancement ratios as high as 3.5 and : 
under | vacuum and ас atmospheric pressure, 
respectively, were obtained. 


The enhancement ratios for finned tubes are 
Ва Боас топат O the tube thermal о For 
he tubes with different thermal conductivity, the 
highest enhancement was obtained for the copper tube, 
while the stainless steel tube had a performance even 
less than the smooth tube under vacuum conditions. 
This is due to_the high wall thermal resistance of the 
stainless steel tube. 
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В. 


RECOMMENDATIONS 


Attach drainage strips on some of the existing tubes 
to investigate the effect of the strips on the 
condensate retention angle and the vapor-side 
heat-transfer coefficient and compare with the already 
existing data. | 


Take data on a tube with a more nearly, parabolic fin 
БАР си and compare them with the existing data of 

ins with triangular, trapezoidal and rectangular 
profiles. 


Take data with different tube diameters to investigate 
its effect on the heat-transfer performance. 


Take data with different. vapor velocities to 
investigate the effect of vapor shear оп the 
heat-transfer performance. 


Take data using different fluids _to study the 


dependence of the heat-transfer coefficient on the 
fluid properties. 
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APPENDIX A 
DERIVATION OF THE HEAT-TRANSFER COEFFICIENT OF FLOODED 
REGION FOR WIRE-WRAPPED TUBES 


Assuming one-dimensional heat conduction for a differen- 
tial element dx, the differential heat transfer is given by 
equation (A.1): 





AT 
dQ = sRtxy " Po dx (А.1) 
where 
Eu c D, 2y(x) 
I R(x) = К; Е | (А.2) 


where y(x) is the vertical distance between the wire surface 
and the horizontal tube surface (see Figure A.1), and is 


given by equation (A.3): 


Ев. - В, - (А.3) 


X is the axial coordinate along the tube. Substituting 
eqation (A.2) for XER(x) into equation (A.l) and integrating 


9 


От Do AT D. T 2y(x) [m Zy(x) (А.4) 
~“ 0 + — 


kr m 





Also heat transfer is given by equation (A.5): 
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D 
Q = h r Do x АТ (A.5) 


Combining equations (A.&) and (A.5), the heat-transfer coef- 


ficient can be expressed as: 





1 Du dx 
= = OT D, Um G) 
kf ky 


Condensate 





Figure А.1 Condensate Film Profile for 
Fully Flooded Wire-Wrapped Tube. 
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APPENDIX B 
LISTING OF RAW DATA 


The following pages contain raw data obtained for tubes 


number 6, 17, 36, 38 and 45 thru 71 under vacuum conditions 


and at atmospheric pressure. 
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APPENDIX C 
UNCERTAINTY ANALYSIS 


Experimentally determined quantities are always associ- 
ated with uncertainties owing to the measuring device accu- 
racy, calibration of the device, and the  operator's 
experience. During this thesis effort, numerical data were 
taken and, together with theoretical formulation, the steam- 
side heat-transfer coefficients were calculated. Since the 
devices used during this experiment to read steam tempera- 
ture, inlet and outlet cooling water temperatures, flowrate 
of the cooling water in finally computing the steam-side 
heat-transfer coefficient, the final result may be distorted 
due to the uncertainty propagation during calculations. In 
cases where the final results show large uncertainties, it 
may be unwise to accept the experimental results. The 
uncertainty on a computation can be determined using the 


following equation proposed by Kline and McClintok [35] 
shown below: 


dR 2 dR 2 eo e ө °... oR 2 YA. 
га ||| a es. + + | Н E.1) 
WR e n) и (55 | = | 


where 


Wp is the uncertainty of the desired dependent 


variable 

X]; X55 ..-., K, are the measured (independent ) 
variables 

Wy. Wo, -+-+-+, М ате the uncertainties in the measured 
variables 
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Using program "UNA7" which is listed at the end of this 
appendix the uncertainties associated with various quanti- 
ties during this investigation were obtained. Also, listed 
in this appendix are some of the selected uncertainty evalu- 
ations. А complete discussion on the uncertainty analysis 


used for this experiment is given by Georgiadis [5]. 
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FILE NAME : UNAZ 
PEVISED February 28. 1986 


CUM обе те Си) 

ОЗЕ 

рата 0. 10086091 , 25727 .94369, - 767345 . 8295, 78025595, 81 
рата -9247486589,6.97688Е +11, -2.66192Е +13, 3.94078Е +14 
READ C(*) 

PRINT 

PRINTER IS 701 

PRINT USING "10Х,""ПАТА FOR THE UNCERTAINTY ANALYSIS:°"" 
PRINT 


BEEP 

INPUT ENTER FILE NANE TIlE 

кш USING "10X,""File Name: "ен ви ељ 
ВЕЕ 


D "ENTER DATA SET NUMBER FOR UNCERTAINTY ANALYSIS",Ids 
Р 

INPUT “ENTER PRESSURE CONDITION (0=V,1=A) “Pre 

Pro Preni 

BEEP 

INPUT "ENTER Ci”,Cı 

ASSIGN @File TO Files 

ENTER *File:Ifg.Inn 

IF ІҒа-0 THEN ENTER Е: Је: ра 

IF Ifq*1 THEN ENTER *File:Dd.Dd.Dd 

FOR I*1 TO Ids 

ENTER GFile:Bvol,Bamp,Vtran,Etp,EC*)o) ,Fm.Tci. Tco,Phg,Puater 
NEXT I 

Ет$ =Е (0) 

IF Prc*1 THEN 


BEE 
ERT USING "10X,""Pressure Condition? Vacuum (11 kPa) s 


PRINT USING TOX. "Pressure Condition: ñtmospheric (01 KPa ae 
END IF 

ЕБЕР 

PRINTER IS 1 

PRINT USING Ax. Select Material Code ГЕ 
PRINT USING "6X.""0 Copper 1 Stainless stes] =~ 
PRINT USING "6X.""2 Aluminum 32 90:10 Cu-Ni""" 
PRINT USING 6х. "А nrFil""" 

INPUT Itt 

IF Itt=0 THEN 

ВЕЕР 

INPUT “SELECT GOSTHIN., I=IHICK)" Tuk 

END IF 

PRINTER IS 701 

IF Itt=0 THEN 

Di=.0127? ! Inside diameter of test tube 
Кс= 385 

Окс=10 

IF Iwt=0 THEN 

0о=.0137 

БЕБЕ 

Во=.01905 ! Outside diameter of test tube 
END IF 
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Ја О 
113778 
1320 
025 
1330 
1585 
1340 
1345 
1350 
1355 
1360 
11855 
1370 
1375 
1380 
1685 
290 
396 
1400 
1405 
1410 
1415 
1420 
1425 
1430 
1435 
1440 
1445 
1450 
1455 
1460 


END IF 

IF Itt=i THEN 

Кс= 16 

Dk c=] 

01=.01245 

Do=.0145 

END IF 

IF Itt=2 THEN 

Kc=167 

Dkc=5 

Do=.0137 

01=.0127 

END IF ; 
Ше * THEN 4 
Кс=45 

Окс=2 

Di=.0127 

Do=.0137 

END IF 

IF Itt=4 THEN 

Кс=385 

Dkc=10 

01=.0156 

Do=.0175 

END IF 

01=.01905 

02=.01587 

IF Itt=4 THEN D2=.01905 
PRINTER IS 701 
Ts=FNTvsv(Emf ) 

PRINT USING "10X.""Steam Temperature = "'. ODNSD, (Пеа 


ЕЛЕК” Тс 


1465 
1470 
1475 
1480 
1485 
1430 
1495 
1500! 
1505 
1510 
1515 
1520 
1525! 
1530 
1585 
1549 
1545 
1550! 
1555 
1560 


PRINT USING "10X.""Hater Fiou Rate (7) eS D oD" Ға 
Dtc:i=.01 

Dtco=.01 

BEEP 

Demt=1.0E-6 
Dts=50R(((C(1)+24C(2)4«Emt+3*C(3)4Emf °2+4*C(4)*Emt °3)*Demf) “2) 
Т-(ТсізТсо)/2 ! FILM TEMPERATURE 

UNCERTAINTY IN THE COOLING WATER 

Drho=.5 ! ERROR IN WATER DENSITY 

Dmt=.0044 ! ERROR IN MASS FLOW RATE 

Rho=FNRho(T) ! WATER DENSITY 

Mf=1.04805E-2+6.80932E-3*Fm ! MASS FLOW RATE OF COOLING WATER 
CORRECT MF FOR THE TEMPERATURE EFFECT 
Mf*Mf*(1.0365-1.966544E-3*Tc1*5.252E-5*Tc1'2)/.395434 
BA1-(PI*Di*^2)/4 ! TUBE INSIDE CROSS SECTION AREA 

Ddi=.000025 

DaisPI*Di*Ddi/2 ! ERROR OF INSIDE TUBE CROSS AREA 

COMPUTE THE WATER VELOCITY 

Vu=Mf/(Rho*#A:1) ! WATER VELOCITY 

PRINT USING "10X,""Water Velocity - Z DD E (m/ 


alae e V 


1565! 
1570 
15 (5 
1580 
1585! 
15-50 
1535! 


CORRECT OUTLET WATER TEMP. FOR THE MIXING CHAMBER EFFECT 
IF Inn=1 OR Іпп=5 THEN TcosTco-.004*Vu'2 

ПЕ ОБЕ ОТНЕМ Teo=ieo-(—. 001384 7001+Vu" 2) 

Te(Tci*Tco)*.5 ! FILM TEMPERATURE 

COMPUTE THE ERROR IN WATER VELOCITY 

DvusVy*SOR(CCGDmf /M£)'2*(Drho/Rho)'2*(Dai1/8B12^2) 
UNCERTAINTY IN THE REYNOLDS NUMBER 


^ 


Е 


1600 Mu-FNMu(CT) ! WATER VISCOSITY 

1605 Dmw=6.E-6 ! ERROR OF WATER VISCOSITY 

1610 Res(Rho*Vu*Di1)/Mu 

1615 DresRe*«SQR((G(Drho/Rho) '2*tDvu/Vu) '2*CDdi/Di1)^2*CDmu/Mu) `2) 

1620! UNCERTAINTY IN THE HEAT TRANSFERRED 

1625 Cpw=FNCpu(T) 

1630 QsMf*(iTco-Tci)*Cpu 

1635 Dcpw=8 

1640 Da=Q*SQR((Dmf/Mf) `°2+((Dtco/(Tco-Tci))) "2+((Dtci/(Tco-Tci)))`2+(DcpmZCpu) 2 
) 


1645! UNCERTAINTY IN THE HEAT FLUX 
1650 Di=.000S ! ERROR IN TUBE LENGHT 
1655 Ddo=.000025 


1660 L=.1333S ! CONDENSING TUBE LENGTH . 

1665 Qp=Q/(PI*Do*L) ! HEAT FLUX 

1670 PRINT USING "10X,""Heat Flux = "" 2.ЗБЕ. = H /m 21: 
"а? р 

1675 PRINT USING "10X.""Tube-metal thermal conduc. = "У ОДО ит (W/m. 
КЕ 


1680 PRINT USING "10X.""Sieder-Tate constant 
1685 Dap-Qp*SQR(C(Da/Q)'2*C(Ddo/Do)^2*(D1/L) *2) 
1690 Lmtd-(Tco-Tci)/LOG((CTs- Tci)/CTs- Tco)) 
1695 UosQp/Lmtd ! OVERALL HEAT TRANSFER COEF. 
1700 AisDts*(Tci-Tco)/((Ts-Tci)*( Ts- Tco) *LOGCCTs- Tc12)/(CTs- Tco)0)) 
1705 Rñ2=Dtei/U(TIs-Tci)=*LOG((Ts-Tci)/(Ts-TIco))) 
1710 A3sDtco/((Ts-Tco)*LOG(CCOTs- Tci1)/CTs- Tco0)) 
1715 Dimtd=tLmtd*SQR(A1°2+A2°2+A3° 2) 
1720 Duo=Uo0*SQR( (Dap/Qp) *2+(Dimtd/Lmtd) *2) 
1725 M=Mu 
1750. ами. 1522346 
1735 Kw=FNKw(T1) 
1740 Ac=0. ! INTERSCEPT FROM SIEDER PROGRAM 
1745 L12.060325 ! LENGTH OF UNFINNED LEFT PART OF TUBE 
1750 L22.034925 ! LENGTH OF UNFINNED RIGHT PART OF TUBE 
1755 Рг=Срш*Миш/Кш 
1760 Мише-ҒММиш(Т) 
s M OF INSIDE HEAT-TRANSFER COEFF. 
Ғ-1. 
1775  His(Kyu/D1)*CCi*Re^ .8*Pr^.333*Cf * Ac) 
1780 Dti-40/(PI*Di*(L*L1*Fel*L2*Fe2) Hi) 
1785 Crc5(Muu/FNMuu(T*Dti2))^.14 
1790 IF ñBSO((GIc-Cif)/G(c)>.01 THEN 
1795_ Ct=(GF+C Oo) 8 5 
1800 GOTO 1775 
1805 END IF 
1810 P1=PI*(Di+D1) 
181S B1=(D1I-D:i)*PI*(Di+D1)*.5 
1820 HBM1=(Hi*P1/(Kc*B1)) .5 
1825 P2=PI*(D:+D2) 
1830 B2=(D2-Di)+*PI*(D1+D2)*.S 
1835  M2-(Hi*P2/(Kc*B2))^.5 
1840 FelisFNTanh(M1*L12)/(OM1*L 1) 
1845 Fe2-FNTanh(M2*L2)/((M2*L2) 
1850 DtcsQ/(PI*Di*(L*L1*Fel*L2*Fe2)*Hi) 
1855 IF ABS((Dte-Dti)/Dte)>.01 THEN 1775 
1860 Dkw=.0010 ! ERROR IN WATER THERMAL CONDUCTIVITY 
1865 Dci=.0005 ! ERROR IN SIEDER-TATE COEFFICIENT 
1870 Dpr=.0S5 ! ERROR IN PRANDTL NUMBER 
1875 Dcf-8.E-6 
1880 A4=.14*Dcf/CF 


че Паре Сл 
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SSS Dhry-Hrs*sSORDODEUZKu) 2 20091701) 2%(.98%Пте/Ке)”2%(,333%һртАВБт) 2%(Пс1/С1) 72: 
AG) 

1890! UNCERTAINTY OF OUTSIDE HEAT-TRANSFER COEFF. 

1895 Rw=Do*LOG(Do/Di)/(2*Kc) ! WALL RESISTANCE 

1900 Ho=1/((1/U0)-(Do*L/(Di*(L+L1*Fel+L2*Fe2)*Hi) )-Rw) 

1905 Drw=Rw*SQR( (Ddo/Do) *2+(Dkc/Ke) *2+(Ddo/(Do#LOG(Do/Di)))°2+(Ddi/(Di1*LOG(Do/D 
my) 2) 

1910 A8521/Uo-Ru-(Do*L/(Di*(L*L T *Fet *L2*Fe2) *H1)) 

1915 A6-Duo/(Uo'2*85) 

1920 A7=Dru/AS 

1925 A8=((Do/(Di*Hi) )*(Dhi/H1i))/AS 

1930 PRINT 

1935 Dho=Ho*SOR(AG*2+A7°2+A8B°2) 

1940! CALCULATE THE X UNCERTAINTY IN Ho 

1945 PrhosDho*100/Ho 

1950! CALCULATE THE % UNCERTAINTY IN REYNOLDS NUMBER 

1955 Prre=Dre*i00/Re 

1960! CALCULATE THE % UNCERTAINTY IN MASS FLOW RATE 

1965 Prmt=Dmt*+100/MF = 

1970! CALCULATE THE XZ UNCERTAINTY IN HEAT TRANSFER 

1975 Prqp=Dqp*100/Qp 

1980! CALCULATE THE % UNCERTAINTY ІМ МТО 

1985 Primtd=Dimtd*100/Lmtd 

1990! CALCULATE THE Z UNCERTAINTY IN Rw 

1995 Рггш= Огш+ 1 00/Кш 

2000! CALCULATE THE % UNCERTAINTY ІМ OVERALL HEAT TRANSFER COEF. 

2005 PruosDuo*100/Uo | 

2010! CALCULATE THE Z UNCERTAINTY IN INSIDE HEAT TRANSFER COEFF. 

2015 Prhi=Dhi*100/Hi 

2029 PRINT 

2025 PRINT USING "10X,""UNCERTAINTY ANARLYSIS:""" 

2030 PRINT 

2035 PRINT USING "10X."'" VARIABLE PERCENT UNCERTAINTY""" 
2040 PRINT 

2045 PRINT USING "10X.""Mass Flou Rate, Md “" Z.2D BID 
2050 PRINT USING "10X.""Reynolds Number. Re п а). сре 
2055 PRINT USING "10X,""Heat Flux, q Ар Рта 
2060 PRINT USING “10X,°“Loge-Mean-Tem Diff, LMTD МА РОА У те 
2065 PRINT USING “10X.""Wall Resistance, Rw DD ED. РЕ 
29429 РЕ НЧ USING ~10X.°"Overall H.T.C., Uo о Оса SPP 
2075 PRINT USING "10X,""Hater-Side H.T.C., Hi "SD: OD UP TE 
2080 PRINT USING "10X,""Steam-Side H.T.C.. Ho ро, сте 
2085 END 

2090 DEF ҒММиш(Т) 

2995 А=247.87(1+133.15) 

2100 Muw=2.4E-5*10°A 

2105 RETURN Huu 

2110 FNEND 

2115 DEF FNTanh(X) 

2120 P=EXP(X) 

2125 O=EXP(-X) 

2130 Tann=(P-Q)/(P+Q) 

=) RETURN Tanh 

2140 FNEND 

2145 DEF FNKw(T!1) 

2150 Киш--.92247»Т1%(2,8395-Т1%(1.800/-Т1%(.52577-.073д4ажТ1))) 

2155 RETURN Kw 

2160 FNEND 

2165 DEF FNMu(T) 

2170 Rñ=247.3/(T+133.15) 


Мш=2.4Е-5+ 10А 
RETURN Nw 
FNEND 

DEF FNRho(T) 


КҺһо"949,52946“Тж(.01268-ТеС524582513Е-5- ТТР ТЕСЕДІ 


RETURN Rho 
FNEND 
DEF FNCpw(T) 


Cow=(4.21120858-T+(2.26826E-3-T#(4.42361E-S+2.71428E-7*T)))*1000 


RETURN Cpw 
FNEND 

DEF FNTvsv(Emf) 
COM CE CEA) 
T=C(0) 

FOR је] 10 7 
T=T+C( I) *Emf I 
NEXT I 

RETURN T 

FNEND 
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DATA FOR THE UNCERTAINTY ANALYSIS: 


ime ЕОБУТ Д5 
Pressure Condition: Deum Cdi Pa 
Steam Temperature - 2 Прати ко 
Мег rios hate (7) = $0.00 
Newer Velocity = 4 e (m/s) 
Heat Flux = оо. ЕЕ ЧУЮ o) 
ве тета: thermal conduc. = 285. 0 (W/m.K) 
Seder lale constant = OBES 
UNCERTAINTY ANALYSIS: 

VARIABLE PERCENT UNCERTAINTY 


S 


то 


Mass Fiow Rate. Ма 
Reynolds Number. Ке 
Heat Flix, а 

LSS Means епш Шу ИЙ 
Wall Resistance, Ru 
ШЕГЕР Н. |, С. . Џо 
EET Side H.T.C.. Hi 
keam- Sae R.T. C., Ho 


cr = Гл =— — -— => 
C 
N 


DATA FORMAE UNCERTAINTY ANALYSIS: 


mee Names Р 05145 
<< поп пе ог: Ол О ЛП ЕР а 
Stes erature = "45, 393 (реа С) 
Натег Flow касе (7) = UU 
Water Velocity 5 1,16 (m/s) 
Heat un - DEUS (H/m 2) 
mee mel thernal conduc. = ОБ. (W/mK) 
eder: Е соп апі = ШЕЕ 
UNCERTAINTY ANALYSIS: 
VARIABLE PERCENT UNCERTAINTY 
es Flon Rate., Md E 
Reynolds Number, Re S E 
Неа? Г] а 210 
LEG eam m Dit f. OMID 0.60 
Wall Resistance. Ru D 
Обе MEI Uo ЕЕ 
ЕКЕ ТЕ moe H.T. C.. Hi n 
ит пе N.T. C.. Ho ЕТ 
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рата FOR THE UNCERTAINTY ANALYSIS: 


File Name: ЕОБА226 | 
Pressure Condition: Atmospheric (101 kPa) 
Steam Temperature = =ч (Dea C) 
Water Flow Rate (4) = 80.00 
Water Velocity = 4,40 m sj 
Неа Ни = 1. ПЙЗӘЕТПЕ PO n 
Tube-metal thermal conduc. = 385.0 (W/m.K) 
Sieder оре ап. = 0.0630 
UNCERTAINTY ANALYSIS: 

VARIABLE PERCENT UNCERTAINTY 
ا اق ع‎ his ОЭ 
Reynolds Number, Re | IS 
Hear оу 089 9 
Log Mean lem Diff, Lae 0.40 
Wall Resistance, Ёш p 
Overall Н [а 
Nater olde mHE C NN EE 
Steam IGE m 6 
DATA FOR THE UNCERTAINTY ANALYSIS: 
File Hame: FO058226 
Pressure Condition: Atmospheric (!01 kPa) 
Steam Temperature =. Чоо (Dea C) 
Water FIO TE (7) = 20.00 
Water Velocity = PIS (m/s) 
Heat Flux = 65.036E*05 (W/m 2) 
Tube -me t3 Mermal Condic. = ЕЭК (W/m.K) 
= тедер=— |а ге Ере рате = 0.0530 
UNCERTAINTY ANALYSIS: 

VARIABLE PERCENT UNCERTAINTY 
Mass Flow Rate. Md 3.00 
Reynolds Number, Re cM 
Heat Flux. д 3.04 
Loa- Mean TED Оена 
Wall Resistance. Rw ЗЕЙ 
Overa Mink Te Б а 
Water-Side H.T.C., Hi RA 
сала” ы I8 
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